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FOREWORD 


This document Is part of a three volume report prepared under NASA Ames 
Contract NAS2-10294, Mathematical Modeling for Vertical Attitude Take-Off and 
Landing (VATOL) Simulation. Volume I: Model Description and Application 

provides background and details of a generic mathematical model for simulation 
of VATOL aircraft concepts. A six-degree-of-freedom off-line (non-piloted) 
digital simulation program incorporating this model was developed and applied 
to the Vought SF-121 VATOL concept. Volume I gives results of this 
application which included development and demonstration of a control system 
for terminal VATOL operations. Volume II: Model Equations and Base Aircraft 

Data gives all the model equations and SF-121 aircraft data in a simulation 
data package format. This volume facilitated the development of a piloted 
VATOL simulation at NASA Ames. Volume III: Users Manual for VATOL Simulation 

Program provides a description of the six-degree-of-freedom off-line digital 
simulation program, instructions for its application, and examples of set-up 
decks and output for several of the SF-121 application runs. 

Project Monitor for NASA Ames was Mr. Gary Kill. The Principal 
Investigator for the Vought Corporation was Robert L. Fortenbaugh. 
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Trim button for fore-aft control stick 
deflection 
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5 


Pedal deflection trimmer 
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Damping ratio of complex root 
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ATT 


Roll attitude feeoback gain (rad/rad or 
rad/sec/rad) 


Roll rate feedback selector gain 
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Roll input proportional gain (rad/unit 
cockpit controller deflection or 
rad/sec/unit cockpit controller deflection) 


Roll input integral gain (rad/ sec/ Pg^ 
or rad/sec 2 /P CMD ) 

Roll control system proportional error 
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Roll control system error integral gain 
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p 
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P CrtD 
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Pc 
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Output of input integrator path (rad or 
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P 


e 


Roll control system error signal (rad or 
rad/sec) 



Limiter on error integrator (units of 
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Output of error integrator path (units of 
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P 
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P TRM 


Alternate name for p c ' when the roll 
attitude feedback signal is Pj N j. (i.e. 
the system controls Pj N j) (rad) 

Trim input from pedal trim or lateral 
stick trim (rad or rad/sec) 


6 R0LL Normalized control output of roll control 

system; roll control input for actuation 
system 

Troll Time constant command filter (sec) 

Logic switch related to CS$ W required to 
insure proper initial ization of roll rate 
integrate, and system trim input. 
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Pitch input integral gain 
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Pitch control system proportional error 
gain (apj-j-^/rad/sec or 6p ITCH/rad) 


K 


q 


e 


II 


Pitch control system error integral gain 
(rad/sec or rad/sec^) 



Pitch rate feedback gain (rad/rad/sec or 
rad/sec/rad/sec) 




Longitudinal stick trim input gains 
(rad/sec/unit input or rad/sec^/unit 
input) 


V 


Pitch attitude feedback selector gain 

Integral of pitch rate feedback selector 
gain 

Input to command filter (rad or rad/s^c) 
Output ot command filter (rad or rad/sec) 


q 


C TRM 


Output of input integrator path (rad or 
rad/sec) 
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q e 


Pitch control system error signal (rad or 
rad/sec) 


q e 


LIM 


Limiter on error integrator (units of 

6 pitch) 


q 


e TRM 


Output of error integrator path (units of 

& pitch) 


q 


I NT, 


q TRM 


Alternate name for q c ' when the pitch 

attitude feedback signal is qj^ (i.e. 
the system controls dj N j) (rad) 

Trim input from longitudinal stick trim 
(rad or rad/sec) 


fipUCH Normalized control output of pitch control 

system; pitch control input for actuation 

system 


TpiJCH Time constant of command filter (sec) 

Logu switch related to CSsw required to 
insure proper initialization of pitch rate 
integrator and system trim input. 
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Lateral acceleration feedback gain 
( 4 y AW /sec 2 ) 
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Yaw input integral gain (rad/sec/R^p or 
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Yaw control system proportional error gain 
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Yaw control system error integral gain 
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O 

Pedal trim input gain (rad/sec /unit 
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(rad/sec 2 /unit input or rad/sec/unit 
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Roll attitude feedback gain to yaw control 
system (rad sec/rad or rad/rad) 



LIST OF SYMBOLS (Continued 


Yaw Control System Symbols (Figure 2-29): 


Yaw attitude feedback selector gain 
Integral of yaw rate feedback selector gain 
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Input to command filter (rad or rad/sec) 
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Yaw control system error signal (rad or 
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Limiter on error integrator (units of 
S YAW) 

Output of error integrator path (units of 

6 yaw) 

Trim input from pedal trim or lateral 
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Alternate name for r c ( when the yaw 

attitude feedback signal is r^y (i.e. 
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Heave control system proportional error 
gain (lb/ft/sec) 
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SUMMARY 


Vouglit Corporation has developed under contract NAS2-10294 for NASA-Ames 
Research Center a mathematical model of a high performance airplane capable of 
vertical attitude takeoff and landing (VATOL). The model is to be applied in 
piloted simulation studies at Ames of the flying qualities of VATOL aircraft 
in terminal area operations. An off-line (non-piloted) digital simulation 
program incorporating this model has been developed to provide trim conditions 
and dynamic check runs for the piloted simulation studies and support dynamic 
analyses of proposed VATOL configuration and flight control concepts. The 
final report for this contract is contained in three volumes. Volume I 
provides development details for the various simulation component models and 
describes the application of the off-line simulation program, VATLAS (Vertical 
Attitude Take-Off and Landing Simulation), to develop a baseline control 
system for the Vought SF— 121 VATOL airplane concept. Volume II contains the 
simulation math model equations and baseline SF— 121 data provided to Ames for 
the piloted simulation studies. The equations are the same as those 
implemented in VATLAS. Volume III is a users manual for VATLAS. It describes 
input data for each program option, contains several illustrative example 
runs, and includes complete program listings. 

The component models for VATLAS were assembled from various sources. 
Several are "off the shelf". For example, the aerodynamics model is a 
modified version of a deterministic, high angle of attack, and large sideslip 
angle model developed at the Naval Air Development Center for application to 
low speed flight of V/STOL airplanes. Likewise, the propulsion system 
dynamics, generalized flight control system, and actuator models were applied 
in the simulation for the recent Lift/Cruise Fan V/STOL Flight Control/Flyirg 
Qualities Program. Other models are new or highly modified versions of 
existing models. Examples of these are the inlet ram effects, Coriolis 
effects and propulsion system-Reaction Control System interactions models. In 
the aircraft equations of motion, the standard Euler rate equations were 
replaced with a direction cosine formulation to avoid the singularity in the 
Euler equations when pitch angle equals 90 degrees. A collection of 
"pseudo-pilot" functions was developed to enable the VATLAS user to impose 
open loop time varying inputs (e.g. steps, doublets) on each cockpit 
controller and/or specify pilot-controlled (closed-loop) stationkeeping and 
transitions. 
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The application of VATLAS to develop a baseline PCS (flight control 
system) for the SF— 121 airplane follows a general FCS design procedure. The 
inputs to this procedure include basic aircraft geometry, mass, inertia, 
lifting surface and fuselage aerodynamics, and propulsion system data; 
operational considerations for the airplane; and a preliminary FCS 
configuration. The operational consideration having most impact on the FCS 
design was that a tilted cockpit (similar to the "nutcracker" concept) would 
be used during hover. This required reversing of the conventional lateral 
stick and pedal roles for hover flight. 

With all inputs defined, a candidate interface between FCS control 
variables and SF— 121 actuators was defined. The acceptability of this 
interface was established by evaluation of six degrees of freedom trims of 
wings level transitions and of hovers in 35 kt winds from any direction. Next 
the FCS gains and control laws were developed by applying linear system 
analysis techniques to the preliminary FCS configuration. Complete design 
analyses were performed at six airspeeds between 0 and 200 kt. Features of 
the FCS thus developed include: 

o Two modes - conventional mode above 60 kt airspeed, hover mode below 60 
kt 

o Three controlled degrees of freedom in conventional mode - roll, pitch, 
yaw; four in hover mode - roll, pitch, yaw, heave 

o Gains are functions of airspeed 

o Attitude stabilization in the rotational degrees of freedom; rate 
stabilization in heave 

o Rate command/attitude hold type system in roll for both modes 

o Rate command/attitude hold type system in pitch conventional mode; 
attitude command in hover mode 

o Automatic turn coordination in yaw conventional mode; attitude command 
in hover mode 

o Rate command in heave in hover mode 
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To evaluate the SF— 121 FCS performance, nonlinear time histories were 
calculated by VATLAS for ten different test cases. The cases Included 
longitudinal stick doublet, lateral stick pulse, and pedal step at 120 kt; 
longitudinal and lateral stick doublets, pedal pulse, and heave rate 
controller doublet at 10 kt; mode switching transient; and pseudo-pilot 
controlled transition and stationkeeping. The response traces are examined 
and discussed in some detail and support the conclusion that the baseline FCS 
provides an adequate starting point for piloted simulation studies and does 
not have to be iterated to correct unanticipated control deficiencies. The 
most distinguishing flying qualities characteristic of the SF— 121 in hover 's 
the non-minimum phase response of the airplane (i.e. it initially accelerates 
in the wrong direction) produced by the use of aft end thrust deflection for 
surge or sway control. This characteristic has the potential for producing 
pilot-induced oscillations when the pilot attempts tight control of surge or 
sway position. 
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1.0 INTRODUCTION 


Aircraft capable of vertical attitude takeoff and landing (VATOL) may 
result from application of certain emerging fighter performance - enhancing 
technologies. These technologies include thrust vectoring, post stall 
controllability, and reliable high authority fly-by-wire control systems. 

VATOL aircraft have many unique handling characteristics, control system and 
display requirements, and pilot-cockpit interface requirements, particularly 
during terminal area operations. A simulation capability is required to 
evaluate and develop these characteristics and requirements prior to a large 
scale commitment of funds for VATOL aircraft development. This capability 
must therefore be available early in the design process as a medium for 
comparing VATOL concepts. Detailed data bases will generally not be available 
for these comparisons. Thus the simulation must be capable of generating 
appliable data sufficiently detailed to emphasize differences and trends among 
several competing concepts. Under contract NAS2-10294 with NASA - Ames 
Research Center, Vought Corporation has developed a generic VATOL simulation 
math model which provides this capability. 

This math model has been implemented at NASA - Ames for piloted simulation 
studies and in an off-line digital prog, named VATLAS (Vertical Attitude 
Takeoff and Landing Simulation). VATLAS provides trim conditions and dynamic 
check runs for the piloted simulation studies and supports nonlinear dynamic 
analyses of proposed VATOL configuration and flight control concepts. 
Development details of the aerodynamics, propulsion system, flight control 
system, inlet ram, and Coriolis effects models and results of a VATLAS 
application to the Vought SF-121 VATOL concept are described in this volume 
which is the first of a three volume final report. Volume II details the 
model equations and b-seline data for the NASA - Ames simulation studies. 
Volume III is a user manual for VATiAS. 

Section 2.0 of this volume details the component models and axis systems 
employed in the complete VATOL model. Interfacing of the models for 
simulation is depicted in figure 1-1. Also described in Section 2.0 is a 
collection of "pseudo-pilot" functions which was developed to enhance the 
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Figure 1-1. VATOL Simulation Math Model General Arrangement 

















utility of VATLAS. These functions enable the VATLAS user to impose open loop 
time varying inputs (e.g. steps, doublets, etc.) on each cockpit controller 
and/or specify pilot controlled (closed loop) stationkeeping and transitions. 

Section 3.0 details the VATLAS application to the Vought SF-121 VATOL 
airplane depicted in figures 1-2 and 1-3. This application is the development 
of a baseline flight control system for the SF-121 beginning with a definition 
of the basic airplane data and proceeding through actuator input specification 
and trim characteristics evalution to specification and evaluation of the 
baseline system gains and control laws. 

Conclusions and recommendations established during the model development 
and application are incorporated in Section 4.0. 
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2.0 VATOL SIMULATION MATH MODEL 


This chapter provides background and details for the six component 
models of the VATOL simulation math model. The component models and 
applicable report sections are as follows: 

o Aerodynamic Forces and Moments - Section 2.2 
o Propulsion System Model Including dynamics, forces and moments, 
and RCS (Reaction Control System) Interactions - Section 2.3. 
c Inlet Ram Forces and Moments - Section 2.4 
o RCS Forces and Moments - Section 2.5 
o Coriolis Forces and Moments - Section 2.6 
o Actuation System - Section 2.7 
o Flight Control System - Section 2.8 

The model discussions are preceded by a description in Section 2.1 of the 
primary axis systems used in the simulation. 

Since VATOL airplanes routinely maneuver at pitch angles approaching 
90 degrees where the standard Euler angle transformation has a singularity, 
the airplane equations of motions were formulated with direction cosines. 

This formulation and other possible alternatives to avoid the singularity are 
discussed in Section 2.9. Section 2.10 details the pseudo-pilot functions 
which were incorporated into VATLAS to simulate pilot control of the aircraft. 
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2.1 Primary Axis Systems 


Four primary axis systems are employed in the VATOL simulation math 
model. These include a North oriented Inertial or earth axis system, the 
vehicle-referenced body axis system, the wind axis system, and the stability 
axis system. The relation between earth axes (Xj, Yj, Zj) and body axes 
(Xg, Yg, Zg) is depicted in figure 2-1. The orientation of the body 
axes with respect to earth axes is determined by an ordered rotation through 
the standard Euler angles - first, yaw (40 around the Z, axis, then pitch 
(©) around the axis, and finally roll ( 4 ) around the Xg axis. Aircraft 
eg position (F ) is measured in earth axes and has components X e , Y e , 
and Z g along the Xj, Yr, and Zj axes respectively. Similarly, 

aircraft eg velocity (r ) Is measured in earth axes and has components 

• • • 

X e , Y e , and Zg. In the body axis system, aircraft velocity has 
components u, v, and w along the X R , y r , and Z R axes respectively. 

X e , Y e , Zg and u, v, w are related by the Euler transformation matrix: 


u 

v 

w 


cos © cos f 
sin 4 sin © cos ^ 
-sinV'cos 4 
cos JPcos 4 sin o 
♦ sin^'sin 4 


cos o sin f' 
sin^sin © sin 4 
♦cosiPcos 4 
sinFcos 4 sin e 
-cos)psin 4 



d 

d 

d 


11 

d 21 

d 31 

12 

d 22 

d 32 

13 

d 23 

d 33 


where the d^'s are direction cosines. 


The body axis system is oriented in the aircraft as follows 
(figure 2-2): The origin is at the aircraft eg. The Xg axis is parallel to 

the fuselage reference line and positive forward. The Yg axis is 
perpendicular to the aircraft plane of symmetry and is positive to the right. 
The Zg axis is in the plane of symmetry, perpendicular to the Xg and Yg 
axes, and positive downward. This orientation remains fixed in the aircraft 
for all time. Locations of aircraft features of interest (e.g. engine inlets 
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and exits, lifting surface centers of pressure) In the body axis system are 
expressed In aircraft coordinates - fuselage station (FS), butt line (BL), and 
waterline (WL). The directions of increasing aircraft coordinates are 
indicated on figure 2-2. 

Unlike the body axis system, the wind and stability systems do not 
maintain fixed orientation to the aircraft. The relations between wind, 
stability, and body axes are depicted in figure 2-3. The X wind axis (Xy 
always points into the relative wind. Note that if angle of attack (a) and 
sideslip angle (b) are both zero, the wind, stability, and body axis systems 
are coincident. Similarly if a » 0 and b * 0, the body and stability axes are 
coincident and, if a ^ 0 and b ■ 0, the stability and wind axes are coincident. 

The VATOL math model uses several axis systems of convenience for 
force and moment producing elements of the aircraft. These systems are 
introduced as required in the model development. Examples include the inlet 
and Coriolis force and moment axis systems. Regardless of the axis system 
adopted for a component model, all forces and moments are resolved eventually 
into body axes. 
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Figure 2-3. Systems of Wind, Stability, and Vehicle Body Axes 
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2.2 Aerodynamics Model 


Operating with the guideline to produce an easily modified model, the 
NASA contract requirement to produce a component i zed aerodynamics model, and 
the fact that VATOL aircraft will experience large angles of attack and 
sideslip; Clark's model (reference (a)) was adopted as a base for the VATLAS 
aerodynamics model. This model “builds" the total airplane aerodynamics from 
fuselage and lifting surface contributions i.e., it is component i zed. The 
model provides continuous aerodynamics functions for all o's and b's 
i.e., -180 £ a < 180 degrees and -90 £ b < 90 degrees. The model is based on 
OATCOM (reference ( b ) ) techniques and thus requires only a modicum of data 
(e.g. lifting surrace geometry, linear lift curve slopes, stall angles, etc) 
to completely model airplane aerodynamics i.e., it is easily modified to 
represent a range of VATOL concepts. In recognition of the limitations of 
OATCOM in predicting aerodynamic characteristics at high a's and fl's, the 
reader is reminded that, in terminal operations, large o's and b's are 
encountered only at low speeds where aerodynamics do not affect significantly 
the aircraft flying qualities. 

Adapting Clark's model to current purposes was guided by a 
self-imposed requirement to obtain reasonable fits of the published Vought 
SF— 121 aerodynamics data in reference (c). Obviously these published data 
could have been loaded directly into the math model, but at the sacrifice of 
overall VATOL model flexibility and timeliness. The strength of Clark's 
approach is its inherent ability to produce consistent aerodynamic estimates 
for aircraft concepts not yet possessing wind tunnel data bases. 

The procedure outlined below was followed in applying and eventually 
modifying Clark's model for VATOL simulation: 

1. Model parameter data were developed for the SF— 121 either from 
the data in reference (c) or from DATCOM. OATCOM was used if 
the data could not be extracted from reference (c). 

2. Clark's model as published was applied to the SF— 1 21 using the 
model parameter data determined in step 1. 
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3. Calculated and published SF— 121 data were compared to indicate 
what model modifications might be necessary to improve its 
predictive capability. 

4. The candidate model modifications were made, evaluated, and 
incorporated as required. 

Equations for the aerodynamics model which evolved from this procedure are 
given in Section 2.0.1 of Volume II of this report. Model parameter data for 
the Si -121 are given in Section 3.3 of Volume II. Results of the published 
and calculated SF— 121 aerodynamics data comparisons and the modifications made 
to Clark's model to produce the VATLAS model are presented and discussed in 
the next section. 

2.2.1 Data Comparison and Aero Model Modifications 

Figures 2-4 through 2-10 compare the published SF-121 data (reference 
(c ) ) with the "initial" and "final" aero models. The SF-121 data are a 
composite of results from various Vought and NASA wind tunnel tests on similar 
configurations. The "initial aero model" results represent Clark's model 
applied to model parameter data developed from references (b) and (c). The 
"final aero model" results represent a modified version of Clark's model 
applied to the initial aero model parameter data base augmented by the 
additional data required by the modifications. 

Comparison with SF-121 data of the results of applying the initial 
model to the SF-121 led to the following observations: 

1. Predicted and SF-121 lift coefficient (figure 2-4) data are in 
good agreement for all a £ 90 degrees. 

2. Predicted and SF-121 drag coefficient (figure 2-5) are in good 
agreement below a ■ 50 degrees. Predicted drag coefficients are 
considerably higher than SF-121 data for a > 50 degrees. 

3. For a < 20 to 30 degrees, the model provides reasonable 
predictions of SF-121 pitch moment coefficient (figure 2-6) 
data. For o > 30 degrees the model predictions and data have 
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Figure 2-5. Comparison of SF-121 Drag Coefficient Published Data With That Calcined by Aero Model 
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Figure 1-6. Compar ;on of SF-121 Pitch Moment Coefficient Published Data With That Calculated by Initial Aero Model 








opposite trends; the model predicts an Increasing (unstable) 
pitch moment as a Increases while the data show a decreasing 
(stable) pitch moment with Increasing a. The apparent cause and 
resolution of this discrepancy are discussed in the next 
paragraph. 

4. Predicted and SF-121 C ] and C Y (figure 2-7) data are 
In reasonable agreement 8f or all 8 a < 90 degrees. 

5. 0 "cted and SF-121 C n (figure 2-7) data are In good 
agreement for a < 20 to 8 25 degrees. For a > 25 to 30 degrees* 
the SF-121 data are considerabley more negative than the model 
predictions. 

6. Below a » 25 degrees, the predicted and SF-121 rudder derivatives 
(figure 2-8) are in reasonable agreement. The prediction 

is reasonable to a * 90 degrees. C nfi and C 
predictions are not good in the a « 36 to 60 
data fall off faster than the predictions. Beyond a * 60 
degrees, the C nfi and Cy 4 predictions are reasonable; 
data and predictions indicate small values of the derivatives. 

7. SF-121 C^ 6 data (figure 2-9) are reasonably predicted for 

a < 30 degrees. At high a, the data indicate aileron reversal 
which is not predicted by the model. 

8. SF-121 C nj data (figure 2-9) indicate a proverse yaw 
tendency f8r a < 15 to 20 degrees which is not predicted by the 
model. The data trend toward adverse yaw with increasing a (for 
a < 35 degrees) is predicted by the model. With the exception of 
more adverse yaw displayed by the data for a » 55 to 70 degrees, 
the model predictions and data have reasonable agreement for a > 
35 degrees. 


degrees range; the 


The result of most concern from ti.„* initial aero model “valuation is 
the poor prediction of pitch moment coefficient for a > 30 degrees. Based on 
the good agreement in lift coefficient data and predictions, it is obvious 
that some, if not all, the pitch moment discrepancies could be resolved by 
allowing an aft shift of the center of pressure (cp) of the aerodynamic 
components with increasing a. Clark's model has no provision for a cp shift 
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Figure 2-7. Comparison of SF— 121 Static Lateral Stability Derivatives Published Data With That Calculated by Aero Model 














Figure 2-9. Comparison of SF-121 Aileron Derivatives Published Data With That Calculated by Aero Model 






with a. Section 4. 1.4. 3 in DATCOM describes a model for the cp shift of 
lifting surfaces. This model was appended to Clark's basic equations. 
Equations were also added to allow the cp of the fuselage to change with a. 
This fuselage cp shift is not covered in DATCOM but was added to compensate 
for residual pitch moment discrep*. ies between model and data after the 
effect of lifting surface c : , shifts are included. In the absence of a data 
base, the fuselage cp shift would not be implemented. Figure 2-10 presents 
the data comparison obtained with the final pitch moment model i.e., Clark's 
model modified to include cp shifts with a for both fuselasge and lifting 
surfaces. The agreement between pitch moment coefficient data and prediction 
is now almost as good as that for the lift coefficient. 

Examination of the difference between drag coefficient data and 
prediction for a > 50 degrees led to a modification in the calculation of high 
a lifting surface drag. Clark's high a drag model fairs lifting surface drag 
coefficient to 1.2 at a * 90 degrees. This limiting value of drag coefficient 
is based on the reference area of the surface. It seems reasonable that, in 
consonance with other calculations in Clark's model the limit drag should be 
based on exposed surface rather than reference area. This modification (i.e. 
multiply 1.2 times the ratio of exposed area to reference area in the high a 
drag equation) was made in Clark's model and produced the much improved high a 
drag coefficient prediction shown in figure 2-5. 

The differences between C ni and Cy 3 predictions and data in 
the a = 30 to 60 degrees range (figure 2-8) are attributed to a gradual, and 
eventually total, loss of vertical tail effectiveness due to masking by the 
wing and fuselage (refer to SF— 121 3-view on figure 1-3). Since Clark's model 
has no provision for this phenomenon, it was incorporated in a table format of 
vertical tail effectiveness as a function of a. The table was quantified 
using Cy fi data which, unlike other side force derivatives, have only a 
vertical tail contribution. For aircraft not possessing a data base, 
reasonable estimates of vertical tail effectiveness can be made by scutinizing 
relative locations of the vertical tail and other aerodynamic contributors of 
the aircraft as a function of a. The much improved predictions of Cy fi 
and ^na with the modified model are shown on figure 2-8. Note that the 
prediction of C^ fi is only slightly affected by this modification. 
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The modification of vertical tail effectiveness also affects the 
sideslip derivatives. The dashed curves for C n and Cj on figure 2-7 
reflect this modification. The Cy dashed curvl Includis this plus 
another modification to be dlscussid below. The prediction is 
Improved at high a with little change In the alreadf good low a prediction. 

The C n prediction is improved only In that, like the data with increasing 
a, the 8 a1rcraft does not reacquire a positive C n once it becomes 
negative. To predict the large C n shown by the*data for a > 30 degrees 
would require modifications to the^model which are highly configuration 
specific and are thus beyond the scope of intended model application. The 
required changes were not made for the SF— 121 application reported herein. 

The other modification reflected in the dashed C y curve on 

figure 2-7 is the addition of a C Y term in the wing contribution 

8 - *» 

l l 

to C y . Calculations of C y with only the vertical tail effectiveness 
modification showed a degraded predlcton capability from that of the basic 
Clark model. This indicated a need for further model refinement. Examination 
of unpublished wing-body C Y data for the SF-121 justified the inclusion 

of a Cy term in the model and allowed its value to be estimated. 

\ 2 
U L 

Section 5. 1.1.1 of DATCOM has a methodology ‘o determine its value in the 
absence of data. The term was incorporated into the aero model and Cy 
was recalculated. As indicated by the dashed Cy curve on figure 2-7 ?he 
model prediction capability has not been enhanced significantly over that of 
the basic Clark model. The modified model is, however, consistent with DATCOM 
methodology and the vertical tail effectiveness change with a. 

Similar to sideslip derivatives, the highly configuration specific 
nature of high a lateral control derivatives of plain flapped elevons (used on 
the SF-121) precludes reliable predictions. Thus the prediction of Clark's 
model for C na and C^ a (figure 2-9) cannot be improved without data. 

d 

The "final aero model" was used to generate all SF-121 aero data for 
the application described in Section 3.0 of this volume. The major 
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discrepancies remaining between data and model predictions are high a C 

and c la * ^ is postulated that the effects of these discrepancies on 

aircraft flying qualities will be minimal because, in SF— 121 terminal 
operations, high o's are encountered only at low airspeeds and low dynamic 
pressures. To prove this contention is beyond the scope of the contract. 

The modifications to the "initial aero model" to produce the "final 
aero model" are reiterated below: 

1. Incorporation of DATCOM method to calculate cp shift of lifting 
surface as a function of o. Also a fuselage cp shift has been 
incorporated to fine tune the pitch moment coefficient prediction 
if supporting data are available. 

2. The o * 90 degrees lifting surface drag coefficient (= 1.2) has 
been referenced to exposed surface area. 

3. A table to incorporate vertical tail effectiveness as a function 
o f a has been made available. 

4* A Cy term has been added in the wing contribution to 
B C 2 

sideforce calculations. 

The equations for the final model are given in Section 2.0.1 of Volume II 
while model parameter data for the SF— 121 are given in Section 3.3 of the same 
volume. 

2. .2 Aerodynamic Rotary Derivatives 

The aero model calculates rotary derivatives. Reference (c) has no 
rotary derivative data for the SF —121 , thus a comparison of model predictions 

and data is not possible. For completeness, a sampling of calculated SF— 121 

* 

rotary derivatives is given in figures 2-11 through 2-13. These were, of 
course, used in the SF —1 21 application described in Section 3.0 of this volume. 


27 




Figure 2-11. SF-121 Pitch Rate Derivatives Calculated by Aero Model 
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Figure 2-12. SF-121 Rol 
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Figure 2-13. SF-i21 Yaw Rate Derivatives Calculated by Aero Model 
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2.3 Propulsion System Model 


The basic propulsion system model for VAT9L simulation has six 
primary distinguishing features: 

o Thrust deflection in two directions via a gimballed nozzle 
arrangement 

o Engine installation angle 

o Nonlinear dynamics - time constant and rate limits vary with thrust 
level 

o Continuously modulated afterburner which can be lit at any throttle 
setting 

o RCS bleed capability 
o One or two engines 

The model is presented as three submodels: The relations for resolving 

individual engine thrusts into body axis forces and moments are developed in 
Section 2.3.1. The engine performance and dynamics model is discussed in 
Section 2.3.2. The modeling of RCS effects on engine performance is discussed 
in Section 2.3.3. 

Data for the SF-121 airplane propulsion system are provided in 
Section 3.4 of Volume II. The SF-121 uses two scaled MFTF-2800-25-1 engines 
whose performance and physical characteristics are described in reference 
(d). The characteristics of these "paper" engines were developed by Vought 
from a Pratt and Whitney parametric cycle analysis computer program and are 
believed to be attainable for a 1995 IOC (Initial Operational Capability). 

The engine dynamics time constant and rate limits data are the same as used in 
the core engine model of the reference (c) study. 

2.3.1 Direct Thrust Forces and Moments Model 

Figure 2-14 depicts the geometry for resolving the thrust of one 
engine into body axis forces and moments. The thrust is assumed to act at the 
center of the nozzle exit face parallel to the nozzzle centerline. To 
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Figure 2-14. Geometry for Resolving Engine Thrust into Body Axis Forces and Moments 



determine the direct thrust forces and moments requires that the aircraft 
coordinates of the thrust application point and the thrust components in the 
body axis system are known. The thrust application point coordinates will be 
develop d first. 

The VATOL simulation model provides for nozzle rotation (i.e. thrust 

deflection) in two directions about the swivel point and for an engine 

installation angle (o^) relative to the Xg axis. The nozzle rotaions 

provide yaw (^y) and pitch (ey) thrust deflections. Since the nozzle 

swivel point remains fixed in the body axis system, it is convenient and 

necessary to determine the aircraft coordinates UFS aBL cw , &WI C1 J of 

oW sw sw 

the thrust applicaton point relative to the swivel point. All three angles - 

V *T. °T - and the nozzle length (1 n q Z ) influence these coordinates. 
Starting with the engine and nozzle centerline coincident and parallel to the 
Xg axis, the following rotat.^n sequence orients the nozzle relative to the 
body axis system (see figure 2-14). First rotate by o y about the Y sw axis 
which is parallel to the Y g axis, then rotate by about the Z Sw 
axis, and finally, rotate by © T about the Y sw axis. The relationship 
describes this sequence and defines aFS sw , ab£ sw , and aWI_ sw is as 
follows: 



The coordinates of the thrust application point are thus defined as: 


FS TAP * FS SW + aFS SW 

BL TAP = BL SW + aBL SW ( 2 * 2 ) 

WL TAP “ WL SW + aWL SW 

To complete the development, the body axis components of thrust must 
be determined. This resolution is essentially the same as that of equation 
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(2.1) except that the sense of y*y rotation is reversed to reflect the sign 
reversal between the direction of increasing WL and that of the Z g axis: 




COS<Jy 0 

sinoy 


COS^y -sin^y 

F 


COSOy 0 

sinoy 


t cor] 

AY T 

- 

0 1 

0 


siny^y cos^y 

0 


0 1 

0 


0 

_ aZ t_ 


-sinoy 0 

COSOy 


0 0 

1 


-sinoy 0 

COSOy 


0 


where 

aXj, aYj, aZy are the body axis components of direct 
thrust 


(2.3) 


Tcqr is gross thrust corrected for RCS interactions as 
described in Section 2.3.3 and flow turning effects 
as described below. 


The equations for the flow turning correction to thrust are as 

follows: 


where 


COR " T AkPL K AFT 

(2.4) 

A “ Apy + 

h ft *ft, h 1 a ft 2 

(2.5) 

ipy « COS - 1 [cos Oy COS^y] 

(2.6) 


' pppl is gross thrust corrected for RCS interactions only 


’FT, 


and K. 


’FT. 


are constants. 


Apy is the geometric thrust turninn angle defined as the angle 
between the engine centerline and nozzle centerline. 


Equations (2.1) through are the basic relations for the direct 

thrust forces and moments model. They are expanded in full in Section 2.0.2 
of Volume II. 



2.3.2 Engine Performance and Dynamics Model 

The engine performance and dynamics model is depicted in figure 2-15. 
Engine performance is represented by the following model parameters: 


1. Fg is the gross thrust level at idle throttle and is a 
fusion of Mach number (M N ). 

2. Fr is the gross thrust level at maximum rpm at minimum 
afterburner setting and is a function of M^. There is an 
assumption here that the thrust at maximum rpm with no after- 
burner is the same as the thrust at maximum rpm at minimum 
afterburner; actually there is a finite but generally negligible 
difference between these two thrust levels. 


3. 

4. 

5. 


F^ax 9 r oss thrust level at maximum rpm at maximum 

afterburner setting and is a function of M^. 
is maximum engine rpm. 

^MAX 

F RpM is fractional engine rpm and is a function of fractional 
non-afterburning chrust level (Tp'. Tp, a model dynamic 
parameter, is the thrust produced by changing spool rpm 


normalized to F r 
• ^MIN 

b. is the inlet mass flow rate when Tp = 1 and is a 

function of M^. inlet mass flow rate is assumed to be a 
function of Tp only; thus, by inference, inlet mass flow rate 
is not a function of afterburner setting. 


7. Ky.j. adjusts the commanded thrust ’evel (T c ) for the reference 
RCS blee -1 ( p ) level. Thus when model output thrust (T Q ) 
is adjusted for actual bleed (BID) level and BID = B^pp, tne 
thrust applied to the aircraft will equal T . The use of this 
parameter assumes that engine controls will be able to compensate 


for known reference RCS bleed. If the controls are not this 


sophisticated then Kg T can be set to 1.0. 


These pare leters are quantified from steady state engine data. Data for the 
SF-121 engines, which are scaled versions of the Mf TF-2800-25-1 engines 
described in reference (d), are provided in Section 3.4 of Volume II. 
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Figure 2-15. Propulsion Syst.mi Dynamics and Performance Model 












The engine dynamics model for non-afterburning thrust (output is 
Tp) indicated on figure 2-15 is nonlinear and was defined by Pratt and 
Whitney for use by Vought in preliminary Type A V/STOL airplane design 
studies. For the current application a model of afterburner thrust dynamics 
(output is t ab) has been placed in parallel with the basic model. The 
afterburner can be lit according to the logic indicated in the figure. With 
the afterburner lit the more rapid afterburner thrust response will enhance 
the s’ower non-afterburning thrust response which is influenced by the engine 
spool inertia. The equations indicated by figure 2-15 are expanded to full 
detail in Section 2.1 of Volume II. 


The engine dynamics model is particularly suited for parametric 
variations to establish propulsion system requirements and control system 
interface. It is also computationally efficient for real-time simulation. 

The application and validity of the non-afterburning thrust model (in essence, 
the spool dynamics model) has been demonstrated in the simulator studies of 
references (e) and (f). Thus, as a baseline, the same values used in 
reference (e) for the nonlinear time constant (Tp^g) and rate limits 
(Tr and T c 


ip ) of the non-af terburning thrust model used in 
reference (e ) M &e been adopted for the analyses presented herein. In lieu of 
definitive data, the time constant of the afterburner thrust model (T ft g) has 
been set at 0.05 sec to simulate the dynamics of fuel delivery to the 
afterburner. Similarly the afterburner thrust rate limits (T 


AB 


MAX 

AB...^ have been se * high ( = * ^0’) to simulate a non-rate limited 


and 


afterburner 


Also T, 


1 £g has been set at 0.10 and T A g 
this establishes that ?be afterburner will be lit whe£ 

T F (I) by 0.10 (l0Zof F ) or when Tp exceeds 1. The 
will remain lit until T.!?} 1 ?) has droppeJ N lc 0.001 (0.1%of F r ) and 


at 0.001; 

exceeds 
burner 


'AB 

is not more than 0.10 larger than Tr(I) 
IN h 


MIN 


2.3-3 RCS Effects on Engine Performance 


The RCS-propulsion system interactions model is depicted on figure 
2-16. Model inputs are Tp ana the output thrust (T Q ) of the engine 
performance and dynamics model. T q i S to be corrected for RCS effects. 
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Figure 2-16. RCS-Propulsion System Interactions Model 










Model outputs are the corrected thrust (T^pj applied to the aircraft and 
the ratio of RCS bleed air currently available to the maximum available 

The maximum RCS bleed air is available when the inlet mass flow rate 
is maximum which occurs when T F = i. As will be shown in Section 2.5, the 
RCS model applies multiplicative factors to scale down the maximum RCS 
capability; K* is one of these factors. The data for K» for the SF-121 
were taken from reference (d). 

Corrected thrust (T^ppjJ is *ormed by multiplying T 0 by a correc- 
tion factor (K'gj) which is based on actual bleed (BLD) and F^ . The 

first step to determine K' BT i S to calculate the bleed available (B avl ) at 
the current uncorrected thrust level. B ftV1 _ i s then compared with RCS bleed 
required ( ®req ) which is generated by the RCS model. If B RE q > B A y|_ 
then BLD is set equal to B ftVL ; thus in this case K' BT will be set by T 0 
and the RCS is saturated. If, on the other hand, B RE q < b AV i ti, en BLD is 
set equal to B RE q; thus K'gy will be set by RCS requirements and excess 

RCS capability proportional to (B AVL - b RE q) is available. Note also that 

if b r£q = B ref (the reference bleed) then K'gy will equal Kgy (figore 
2-16) and 

^APPL will equal commanded thrust (T c ). The data to define 
K gj for the SF-121 were also taken from reference (d). 

The relations and logic indicated by figure 2-16 are expanded to full 
detail in Section 2.1 of Volume II. 
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2.4 Inlet Ram Forces and Moments Model 


The development of the inlet ram forces and moments model is 
presented in two steps. First the relatively simple, standard model is 
given. This model is phenomenologically correct but lacks experimental and/or 
theoretical guidance to define the ram force application point. This aspect 
of the model s normally left to user whim. A recent Vought contracted effort 
(reference (g ) ) , undertaken to remedy this situation, resulted in a series of 
inlet force and moment design charts intended for preliminary design 
application. Tne standard model was expanded and modified to incorporate 
these results. This is the second step of the model development. 

2.4.1 The Standard Model 

Figure 2-17 depicts the geometry of the standard ram forces and 
moments model in the aircraft plane of symmetry. The forces and moments have 
contributions from two sources: 


1. The forces and moments imposed on the inlet by the captured stream 
of air which enters it. 

2. The additional forces and moments imposed on the inlet when, 
because of aircraft rotation, it moves relative to the captured 
stream of air. 


The model of 
but opposite 
point (AP). 


the f.rst contribution assumes a ram force fF RAM ) parallel 
to the relative airspeed vector (V A ) acting at an application 
The ram force is defined as 


F RAM q * ~ V A 

where m^ is the inlet mass flow rate determined by the propulsion system 
model and V A has the components u^, v^, and w^<j along the Xg, 

Yg, and Zg axes, respectively. Thus the body axis components of F R ^ M 
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Figure 2-17. Geometry of Standard Inlet Ram Forces and Moments Model 


AX R q * “ m I U AS 
aY R q * -®I V AS 


(2.4.1) 


aZ R q - m I W AS 

Anticipating the modifications to the standard model to incorporate 

the recent Vought results, an alternate formulation of aX R t A y R , and 

aZ R is now introduced. The angles A-j-^ rn and Byyp^ are required 

for°this purpose (figure 2-18). Aj URN is defined as the geometric inlet 

flow turning angle and is equal to the°angle V* A makes with the inlet 

centerline. B^ RN is defined as the angle between the V A -engine 

centerline plane and the aircraft plane of symmetry. Using Ayy RN , 

B TURN* and the ram f° rce magnitude (i.e. |F RAM | = my I V A | * ° 

• 0 * 

)» equations (2.4.1) can be recast as 


aX R q a -"i 'U cos a turn q 

aY R o - -“I sin A TURN 0 sin ^TURN < 2 - 4 * 2 ) 

aZr q - -nij V* sin A TURN ^cos &y URN 

where 

A TcRN 0 “ cos ~ 1 ^AS^ 3 

b turn * tan_1 ^ v as /w as3 


As noted above, guidance on where to place AP is lacking. It is 
generally agreed that it should lie on the centerline of the engine. There is 
no agreement, however, in where it should lie relative to the inlet face. It 
has been placed by various researchers from X AppL . 0. upstream to X AppL - 
Djfl and beyond. Since information now exists to quantify X AppL (reference 
(g)), the derivation of the standard mooel continues with X AppL 
unspecified. The moments about the aircraft eg due to F RAM are 

M RAM q - l r AP x ~RAM q ) 
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( 2 . 4 . 3 ) 


l ram 0 


aZ r 0 y in - aY r 0 z in 

H 

- 

aX r 0 z in “ aZ r 0 ( x in + x appl) 

>„ 


aY R 0 ( X IN + X APPL^ " aX R 0 Y IN 


where M|^ m i S the moment vector of 

r A p ?s the vector location of the Xp relative to the eg. 

X IN* y in» Zin are the coordinates of the inlet in aircraft 
body axes. 


l RAM . M ram , N ram are the components of fl RAM 
alon§ the a^rcraft^ody axes. 0 


The second contribution to the Inlet ram forces originates in the 
velocity, produced by aircraft rotation, of the Inlet relative to the Incoming 
stream of air. The defining relation is 


• F RAMj ■ -^i (“a x r IN^ 


~aX B “ 
K I 


“***I (^ Z IN ■ rY IN^ 

aY r 

R I 

m 

-«l( rX IN - P Z IN) 

A 

1 

“^1 ( pY IN ■ qX IN^ 


where 


aX r , aY r , aZ R are 
p, d)» r are the x body 
components of 


the body axis components of F RAM 

axis roll, pitch, and yaw rates ind are also the 


^RAM. acts at the inlet face, 
defined as 


The moments due to 



are thus 
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M RAMj = ^IN X F RAMj ) 



where 


L RAI.j» m RAMj» n RAMj are the boc| y axis components of M RAM ^ 


Combining the two contributions represented by equations 2.4.2 
through 2.4.5 gives the standard model form for the body axis forces and 
moments produced by inlet ram effects: 


X RAM = aX R o + aX Rj = cosA TURN 0 +qZ IN _rY IN^ 

Y RAM = aY R o + = - n, l( V oo s1nA TURN Q sinB TURN +rX IN -pZ lN^ 

Z RAM = aZ R q + aZ Rj = “‘M'to sinA TURN o cosB TURN +pY IN _qX IN^ 

L RAM " L RAM q + L RAMj = y IN Z RAM " Z IN Y RAM (2.4.6) 

M RAM = M RAM q + M RAMj = Z IN Y RAM " X IN Z RAM + X APPL m I ^ sinA TURN 0 cosB TURN 

N RAM = N RAM q + N RAMj * X IN V RAM " X APPL m I sinA TURN 0 sin B TURN “ X RAM Y IN 

Recall that the standard model makes the following assumptions about 

the ram force (Fram 0 ) ; 
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1* I p ram I 1S equal to mj V^ (where ^o«|v A J) for all values 


of a: 


TURN, 


2* Pram acts P^ral lei and opposite to V A 
0 


Vought's recent efforts (reference (g) ) to apply sophisticated inlet modeling 
techniques to determine the ram force application point determined that not 
only the application point but also the actual inlet flow turning angle and 
magnitude of the ram force vary with A JURN , ratio of ambient airspeed to 
Inlet air speed (V*/Vj), and Inlet geometry. The next section will 
describe the modification and expansion of the standard model to incorporate 
these effects. 


2.4.2 Expansion and Modification of the Standard Model 

A brief synopsis of the work performed by Vought in reference (g) is 
required as background for the development in this section: Under contract to 

Naval Air Development Center, Vought has developed a computerized prediction 
method for propulsive induced forces and moments in transition and STOL 
flight. This method is baseo on Vought's V/STOL Aircraft Propulsive Effects 
(VAPE) program. One of the 7APE options provides for the calculation of inlet 
forces and moments. This option combines a highly modified Stockman inlet 
flow model with a program which integrates the inlet pressures determined by 
the inlet model to produce inlet forces and moments. Arbitrary axisymmetric 
inlet geometries are accepted by the pronram. The bulk of the program 
applications in reference (g) is to the NASA QCSEE GE2 inlet. Being a 
subsonic inlet, it has relatively thick lips. To determine the effects of 
thinner supersonic lips, a limited program application was made to a thin 
lipped configuration. The force and moment differences were minor at low 
values of \U/Vj (< 0.3 to 0.4) and negligible at higher values of 
V^/Vj. Model results designated for preliminary design application can oe 
presented as the difference between the geometric and effective flow turning 

angles (^^1, the ram effectiveness factor (R^), and normalized (to 
equivalent inlet diameter) ram moment arm ( x appl^in^ as functions of 
A TURN and W V I (figures 2-19 to 2-21). 
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Figure 2-20. Effect of Velocity Ratio and Geometric F’ow Turning Angie on the Difference 
Between Effective and ^ ornttric Flow Turning Angle 
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Figure 2-21. Effect of Velocity Ratio and Geometric Flow Turning Angle on Ram 
Effectiveness Factor (R M ) 
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There is no provision in the Vought VAPE program for rectangular 
inlets, thus data equivalent to that of figures 2-19 to 2-21 cannot be 
generated for rectangular inlets. For application of the VATOL simulation 
model, it is assumed that the data for axisymmetric inlets can be applied to 
rectangular inlets by replacing them with equivalent circular inlets. The 
equivalent circular inlets have the same area as the rectangular inlets; the 
inlet diameter (D^) specified in the simulation model is the diameter of 
the equivalent circular inlet. The validity cf this approximation cannot be 
readily assessed. It 's believed tha* the inlet forces and moments calculated 
by the simulation model will provide reasonable estimates of actual ram 
effects and will certainly be adequate for comparing inlet ram effects on 
various VATOL aircraft concepts. 


Introducing the data and terminology of figures 2-19 to 2-21 into the 
standard model (equation (2.4.6)) gives the following set of equations: 


A TURN 0 * cos 1 ( u as /v «o ) 

(2.4.7) 

aA TURN - f l RAM tA TURN 0 * t'U/Vi)] 

(2 4.8) 

R M - f 2 RAM £ A TURV ^ /V I^ 

(2.4.9) 

X APPL /0 IN = f 3RAM^ A TURN 0 » (V V I^ 

(2.4.10) 

F RAM = m lW*M 

(2.4.11) 

A TURN = A TjRN + aA TURN 

0 

(2.4.12) 

X RAM * ~ f RAM cos A TURN " *"I ^ Z I N -r Y In) 

(2.4.13) 

Y RAM = _F RAM sin A TURN S1n E i URN ^ rX IN _pZ iN^ 

(2.4.14) 

Z RAM = ~ F RAM sin A TURN cos B TURN " m I IN _qX lN^ 

(2.4.15) 

L RAM = Z RAM Y IN " 'RAM 

(2.4.16) 
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M RAM = X RAM Z IN _Z RAM X IN ~( X APPL /D IN) D IN F RAM sin A TURN cosB TURN 


(2.4.17) 


N RAM * Y RAM X IN ~ X RAM Y IN + ( X APPL /D In) D IN F RAM sin A TURN s1n B T'JRN (2-4.18) 

where Ay^ is the effective inlet flow turning angle 
F RAM t ^ ie effective magnitude of the ram force 

\ 

Equations (2.4.7) through 2.4.18) are the basic relations of the 
VATOL simulation inlet ram forces and moments model. Jhe functional relations 
indicated by equations (2.4.8) through (2.4.10) are the data of figures 2-19 
to 2-21 extrapolated as indicated to (W v j) = 0 ana A TURN = 0 and 180 
degrees. The detailed equations actually programmed in th? VATOL math model 
are given in Section 2.0.3 of Volume II. These Volume II equations include 
the effect of engine tilt angle (a^) which was not introduced here for 
reasons of brevity and clarity. 
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2.5 Reaction Control System Forces and Moments 


The features of the RCS forces and moments model are as follows: 

o Up to ten jets locatable anywhere in the airplane and at any 
angle relative to aircraft axes. 

o Jets can be specified as demand or continuous bleed. Oemand 
bleed jets can be further specified to demand more bleed than 
the RCS reference bleed. 

o Jets thrust in one direction only. 

o The bleed required by the RCS is monitored and limited, if 
necessary, by the bleed available from the engines. 

o Continuous jet forces follow RCS actuator outputs with no lag. 
Demand jet forces are lagged relative to actuator outputs to 
simulate bleed flow dynamics. 

The model is presented in three sections: the procedure for locating and 

specifying the type of jet is described in Section 2.5.1, calculation of 
individual jet forces is developed in Section 2.5.2, resolution of jet forces 
into boay axis forces ana moments is developed in Section 2.5.3. Data for the 
SF— 121 RCS are presented in Section 3.6 of Volume II while detailed equations 
are given in Section 2.0.4 of the same volume. 

The overall approach to RCS modeling is to calculate an ideal force 
demanded at each jet based on an ideal maximum. This maximum is based on the 
maximum inlet mass flow rate of the engines. The ideal force is then 
corrected for actual inlet mass flow rate, bleed available for the RCS, and 
demand for thrust at the other jets. This modeling approach has evolved 
primarily from the author's experience with the math model and background data 
of the VAK-191B airplane (reference ( h ) ) . The procedure for locating jets is 
based on reference (a). 

2.S.1 RCS Jet Parameters 

The parameters involved in locating or orienting RCS jets are 
depicted in figure 2-22. The application point (i.e. location) of a jet force 
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Figure 2 - 22 . Geometry for Locating and Orienting RCS Jets 



is defined in airplane coordinates. The orientation is defined by the ordered 
rotation required to align a jet force initially along the positive Xg axis 
with its final installed force direction. The ordered rotation proceeds as 
follows: first, rotate around the Z g axis by the jet yaw angle (f JET ), 

then, rotate around the rotated Y axis by the jet pitch angle Five 

parameters are therefore associated with the location and orientation of each 
jet: 


o Fuselage Station (FSjgy), Butt Line ( BL ) , and Waterline 
(WLjet) of the force application point. 

o Jet yaw and pitch angles (^j E j and ©jet) 

Three more parameters are required to completely specify each jet. 
These are 5R5, the demand parameter, BLDM, the bleed more parameter, and 

^RCS » the maximum force of the jet. A demand bleed jet is designated 
bysHt ing DMD * 1.; for a continuous bleed jet, DMD » 0. If the jet is 
dem i bleed, then 8 lT5m must be specified; if 6 l DM = 1. the RCS bleed can be 
inct eased beyond the reference level if the additional bleeo capability is 
required and available; if Bl£)M « 0., RCS bleed cannot exceed the reference 
level. (kfiM can be set to 0.0 for continuous bleed jets). The SF-121 has 
t w c RCS demand bleed jets for roll control and has a 3.5^ reference bleed 
level (which is the bleed required for non-RCS uses). Therefore DMD ■ BLDM - 
1. for both SF-121 RCS jets. 

The maximum force of the jet is the thrust produced if the jet area 
is maximum an' bleed required is less than bleed available, the engine is 
operating with maximum inlet mass flow rate, and there is no thrust loss due 
to other operating jets. For the SF-121, F RCS 0 f each jet was set at 
1500 lb. M * 

2.5.2 Calculation of Individual RCS Jet Forces 

Figure 2-23 outlines the calculation of individual jet forces. The 
inputs to the process are the individual normalized jet areas ( 6 rcs ( ^ ) ) • 

T.vse are processed by the standard demand bleed or continuous bleed force 
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(from RCS - Propulsion 







shaping functions to produce commanded forces. The demand bleed shaping 
function Is as follows: 


fRCS^D «RCSU> 


f rcs(^ 


s RCS( I ) > °* 
i RCS^ I ^ < °* 


The continuous bleed shaping function is as follows: 

f rcs( ! ) * f rcs mx ( 1 ) f 1 + 6 rcs( ! ^ 

Continuous bleed jets are assumed to act in pairs. The total jet flow area is 
constant for these jet pairs; differential control forces are produced by 
increasing one jet area at the same rate as the other is decreased. Thus when 
«RCs(I) « 0, the jet areas are equal and both jets command F RC $ ( I } . 

When 6 R qs(I) « 1.0, the one jet area is double its value at fi^fl) * 0» 

and commands 2 (i) while the other jet area is zero and commands 

* MX 

zero force. 


The demand jet commanded force takes three paths (figure 2-23): one 
path determines whether more bleed should be allowed to handle the additional 
RCS force requirement. This additional bleed, if any, is reflected in the 


variable aF 


RCS 


force summatio8 M ?F R ££ ). when all jets have been processed, F R £$ 
represents the total eiut flow area of the RCS. The third path adju^s the 
commanded jet force for various thrust-loss contributors (such as available 
bleed less than required and less than maximum inlet mass flow rate) before 
applying it to the aircraft. The continuous jet commanded force is also 

combined into F R ^ and adjusted for thrust losses before being apolied 

, .. , SUM 

to the airplane. 


The second path combines the force into the commanded 


Bleed required by the RCS (B r ^q) is generated from the sum of the 
RCS force required at the reference bleed level (F R ^ ) and 


aF 


MX. 

0 

Rrs • Note that if the RCS is totally demand bleed, F R rc »0. 
DMD 
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Br£q is compared with bleed available ( ) to establish the actual RCS 
bleed level (BLD). Bleed available is generated by the propulsion system 
model. The relation between RCS force and engine bleed denoted by 
is a function of the engine installation and bleed capability. The variable 


RCS represents the total RCS force available and is less than or equal 
AL 


to F 


RCS, 


KJf, is generated by the propulsion system model and is the 


ratio factual inlet mass flow rate to the maximum available. Thus the 
relation for adjusting the individual jet command forces to the force applied 
to the aircraft can be shown to be: 


F RCS (; ( i ) = Frcs^) K m F RCS AL /F RCS SUM 

This relation assumes each jet is affected proportionally the same by inlet 
mass flow rates and bleed capability. 


F RCS (0 is processed by the RCS force dynamics to incorporate 
the effects of demand or continuous bleeding on the force applied to the 
aircraft. Continuous bleed jets are assumed to have no lag and (i) 
is applied immediately to the aircraft. Demand bleed jets are assumed to have 
first order lag dynamics to represent the delay between jet area changes and 
the appropriate force changes. Thus, for demand bleed jets, F^ (i) i S 
sent through a first order lag before being applied to the aircraft. 


2.5.3 Resolution of RCS Forces into Body Axis Forces and Moments 


Body axis RCS forces are the summation over the number of jets of the 

body axis components of the individual jet forces ( f rCS (U): 

A 


X RCS = 


n JET 

T F RCS cos(fj ET (I))cos(e JET (l)) = 
1=1 A 



aXrcsU) 


n JET 

Y RCS = F RCS.( i ) sin(^] ET (I))cos(oj ET (I)) * 

1=1 


n JET 

27 aY rcs^) 
1=1 
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n JET "JET 

Z RCS - ^ " f rcs( 1 ) sin(o Jet (I)) - Z aZ rc$ (I) 

1-1 A 1-1 

where y rc $, Z RC $ are the total body axis RCS forces 

aX RCS(^» aY RCS( i )» aZrcs(I) are the body axis RCS forces due to jet I 
n JET is the number oY RCS jets. 

Body axis moments are the summation of the moments produced by each jet 


n JET 

L RCS “ 5 t aZrcs ( 1 ) Y jet ( I ) -aY rcs ( 1 ) Z JET ^ ^ 

1-1 


n JET 

m rcs ■ 21 ^ aX rcs( i ) z jet( i ) _aZ rcs( i ) x jet( 1 ^ 

1-1 


n JET 

N RCS * 21 ^Y RCS (I)Xj ET (I)-aX RCS (I)Yj ET (I)] 

1-1 


where l RE <j, M R q^, N R q^ are the total body axis RCS moments 

Xj ET {I), Yj E t(I)» Zjet(D are the boc *y axis com P° n ents of 

the position vector between the aircraft eg and the application point 
of jet I. 
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2.6 Coriolis Forces and Moments 


Forces and moments are generated when the ai craft rotates and 
produces a change in the direction of the mass flow passing through the 
propulsion system. These forces and moments are due to the Coriolis 
acceleration imposed on the aircraft - 'a cor *Cm(w A x 7) where T A is the 
aircraft rotation vector and mV is proportional tp the mass flow rate. 

Although Coriolis forces and moments are not large in magnitude and are 
generally not influential contributors to aircraft dynamics, a model of the 
effect was incorporated for completeness in the WATOL simulation. 

The model is derived from the straight through propulsion system 
representaion shown in figure 2-24. This is an assumption in that the flow 
through actual propulsion systems may have several turns to traverse before 
exiting. The model also assumes that the amount of mass within the propulsion 
system is proportional to the inlet mass flow rate, thus the effect of fuel 
mass flow added within the propulsion system is neglected. This is reasonable 
since fuel mass flow rarely exceeds 10%of the inlet mass flow and is 
typically 5%or less. 


The development of the Coriolis force and moment model equations is 
based on reference (i) and proceeds from the defin'tions of Coriolis moment 
about the eg ( m cqr ) and force (F^qr) vectors: 


~C0R = " 2 |(*A x R ) pA d $ t 

■^ducH ) 


M, 


'COR = -2 f (R X (7 a X R)) pA d^ 

“^ducH) 


where 


( 2 . 6 . 1 ) 


( 2 . 6 . 2 ) 


is the integration variable which moves along the 
propulsion system centerl ine (figure 2-24) i.e. along 
the Xq axis. 
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1 DUCT is 1en 9 th along the centerline from Inlet face to 
exit plane with no nozzle deflection. 

pAR is the mass flow rate through the propulsion system 
R is the vector from the aircraft eg to pointy 

In terms of propulsion system and aircraft variables, the vector quantities in 
equations (2.6.1) and (2.6.2) can be written as 

pAR = -mjT + 0 J + 0 k (2.6.3) 

= (p cos Oy -r sin <jy) i + qj + (r cos + p sin a y )k (2.6.4) 

R = + Yj J -a T (2.6.5) 

where T, J, k are unit vectors along the propulsion system Xq, Yq, Zq 
axes. 

Suostituting equations (2.6.3) through (2.6.5) into equations (2.6.1) and 
2.6.2) gives: 

1 

F C0R = 2 ™I j [( r cos oy + p sin o y )j + q T] dCp (2.6.6) 

" ^OUcH* 

1 

M C0R = ~^\ \ t Y I q ' a ( r cos a y + P sin °y)jl - q J 

- ( 1 ouCT - ^ ) ~ (rcosay + p sinoy)k d^ (2.6.7) 

The Coriolis fore, and moment equations whief resul : when equations (2.6.6) 
and (2.6.7) are integrated are detailed in Secton 2.0.5 of Volume II. 
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c,i ctuation System Model 


The VATOL simulation models for the actuation and flight control 
systems are perceived to in*’rface as shown in figure 2-25. The purpose of 
the flight control system (FCS) model is to combine cockpit controller and 
nation sensor inputs through various control laws to generate control commands 
for the aircraft degrees of freedom. These control commands direct the 
application of body axis control moments and forces to the airplane. The 
act. Hi *n system .-'del converts the control commands into commands for the 
actuators of the control force and moment generators and provides dynamic 
models ijr the actuators. The propulsion system is shown on figure 2-25 
because one of its functions 's to serve as the actuator for thrust commands. 

The four FCS submodels - roll, pitch, yaw, and heave control systems - 
are justified and presented in Section 2.8. The actuation system actuator 
input submodel is discussed in Section 2.7.1 while the actuator dynamics sub- 
model is discussed in Section 2.7.2. 

2.7.1 Actuator Input Model 

The VATOL simulation model has remained generic to this point: 

Insofar as a proposed configuration has no more than four lifting surfaces 
(i.e. left and right wing halves, horizontal stabilizing surface, and vertical 
stabilizing surface), one fuselage, ten RCS jets, and two jet engines with the 
capability to deflect thrust in two directions, it can be simulated without 
changing the force and moment models for any arrangement of these components. 
The futility of maintaining generality in the actuator input model and the 
need for this model to be more configuration specific become apparent in the 
following discussion: The VATOL force and moment models accept inputs f^om 

twenty eight control force and moment generators. As a minimum, the actuator 
input model, to maintain generality, must provide paths from each FCS control 
command to each actuator. Generally the gains in these paths must be 
programmed as functions of aircraft state or control variables to accommodate 
controls blending and/or nonlinear gearings for optimum flying qualities 
during transition. Crossfeeds between actuator inputs are also generally 
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Figure 2-25. Interface of Flight Control System and Actuation System Models 









required to minimize control couplings. An actuator input model which 
provides a priori for all these contingencies would be so large as to preclude 
real time manned simulation. In a practical application, only a small portion 
of the general model would ever be required and specification of that portion 
would be guided by a control system analysis of the simulated aircraft 
concept. Thus, the actuator input model is constrained to be more 
configuration specific than any of the other simulation component models and 
must be developed by off line analyses. 

The development of the SF— 121 actuator input model is described in 
Section 3.5 of this volume. The model is depicted in figure 3-2 and detailed 
by the equations of Section 2.2.1 of Volume II. 

2.7.2 Generic Actuator Dynamics Model 

With the exception of propulsion system thrust dynamics, the dynamics 
of all actuators in the VATOL simulation model are represented by the position 
and rate limited first order model shown in figure 2-26. The position and 
rate limits and time constants for the SF— 121 actuators are specified in 
Section 3.7.2 of Volume II. The model equations are detailed in Section 2.2.2 
of Volume II. 
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Figure 2-26. Generic Actuator Model 





2.8 Flight Control System Model 


As Indicated by figure 2-25, the FCS model has four submodels. These 
are the roll, pitch, yaw and heave* control system models and correspond to 
the four degrees of freedom which normally require some level of automatic 
control or augmentation In VATOL airplanes. The roll, pitch, and yaw degrees 
of freedom are augmented In all flight regimes while heave is augmented only 
In the low speed (hover) regime. Surge and sway control at all speeds and 
high speed heave control are not modeled because these are totally manual 
modes utilizing throttle commands and/or aircraft attitude as the means to 
generate control forces. Generic roll, pitch, yaw, and heave control systems 
are depictec on figures 2-27 through 2-30 and are discussed In Sections 2.8.1 
through 2.8.4 

The FCS model receives inputs (figure 2-25) from three conventional 
cockpit controllers - pedals, control stick with right-left and fore-aft 
degrees of freedom, and manual throttles - plus one additional controller for 
heave. The pedals and both control stick degrees of freedom can be trimned. 
Motion variables assumed available for FCS control law formulation include 
body axis roll, pitch, and yaw rates (p, q, r), stability axis yaw and roll 
rates (r $ and p s ), airspeed (V A ), Euler angles (©, 6,y*). aircraft 
lateral and normal accelerations at the eg (n and n ), and heave 

* r n 


*As used herein, heave, surge, and sway degrees of freedom are referenced to 
earth axes; heave is motion along the Zj axis, surge is motion along the 
Xj axis, and sway is along the Yj axis. 
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The FCS submodels are similar in several respects: All feature the 

option to switch control laws and cockpit controller function as a function of 
the control system switch (CS SW ) value. The value of CS SW is determined 
b y CS$ W = 0. for high speeds, = 1. for hover. Most of the system 
gains are programmable as functions of V A . The forward loops can be set up 
as proportional or proportional -pi us - integral controllers. The cockpit 
controller commands can be input through pure or shaped (by a lag filter) 
gains (for rate or attitude command systems) or through proportional -plus- 
integral arrangements (for rate command-attitude hold or acceleration 
command-rate hold systems). 

Most of the feedbacks provided for the various FCS submodels are 
fairly conventional: For example, body and stability axis roll rates (p and 

P s ) and Euler roll angle (d) are available in the roll control system; body 
and stability axis yaw rates (r and r $ ) and Euler roll and yaw angles (d 
andf') are available in the yaw control system; body axis pitch rate (q) and 
Euler pitch angle (©) are available in the pitch control system; heave rate 
(Z g ) is used in the heave control system. Several unconventional feedbacks 
have also been provided: For example, to avoid the high sensitivity of d 

and V* to aircraft rotations when © > 80 to 85 degrees and provide a feedback 
proportional to attitude, the integrals of roll and yaw rates (Pj N j and 

r INT^ have been mac * e available in the roll and yaw control systems. Also to 
avoid problems with the range of © (-90 < © < 90 degrees) and provide a 
feedback proportional to pitch attitude, the integral of pitch rate (dj N j) 
has been made available in the pitch control system. 

Two areas of concern in VATOL terminal operations which will be 
studied extensively via manned simulation are cockpit controller function 
switching and control mode switching. Cockpit controller function switching 
is required to avoid pilot confusion with VATOL concepts which rotate the 
cockpit as the aircraft approaches hover. The requirement for control mode 
switching should be obvious - the aircraft flies differently in aerodynamic 
force - supported flight than it does in thrust - supported flight. As 
specified in the FCS model, both cockpit controller function switching and 
control mode switching are controlled by the control system switch (CS^). 
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It is asumed that the function and mode switchings occur simultaneously and 
that CSj W changes from 0 to 1 or vice versa In one sample period or one 
computer time frame. Logic Is provided for reinitializing any Input or 
feedback Integrators when CS^y changes value. This helps to avoid large 
switching transients, particularly when switching between Euler angle and 
Integrated rate feedbacks. This switching scenario Is rudimentary and should 
provide a “worst case" starting point. Many switching functions will be 
explored In evolving VATOL design guidelines; each of these will require 
reprogramming of the simulation model. 

2.8.1 Roll Control System Model 


The roll control system model is depicted on figu-e 2-27. Table 2-1 
gives the parameter values required to implement several common combinations 
of roll control system types and forward loop controllers. Feedback variable 
selection is controlled by the three selector gains; K p for roll rate, 
for d, and K/ for p j N j; which are functions of CS$ W and 
generally assume values of 0. or 1.. Those gains which can assume any value 
are indicated in Table 2-1 as being 4 0. The exact values of these gains must 
be established by off-line control system analyses. 


Figure 2-27 represents the roll control system for a VATOL airplane 
in which the cockpit rotates as hover is approaches, cockpit control of roll 
is switched from lateral stick to pedals and roll trim function switches from 
lateral stick trim controller to pedal trim controller. This feature can be 
eliminated from the model by removing the CS$y dependence from the cockpit 
controller inputs to the system. 


Equations for the roll control system model are given in Section 
2.3.2 of Volume II. Parameter values for the SF-121 roll control system are 
given in Section 3.8.2 of Volume II. These parameter values were established 
by the control system analyses described in Section 3.7 of this volume. 
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cokmwuH Figure 2-27. Generic Roll Control System Model 
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2.8.2 Pitch Control System Model 


The pitch control system model is depicted on figure 2-28. Table 2-2 
gives the parameter values required to implement several common combinations 
of pitch control system types and forward loop controllers. Feedback variable 
selection is controlled by the two selector gains; K^ q f 0 r q IN j and K e 
for e; which are functions of CS $W and generally assume values of 0. or 1.. 
Those gains which can assume any value are indicated in Table 2-2 as being £ 0. 
The exact values of these gains must be established by off-line control system 
analyses. 


Equations for the pitch control system model are given in Section 

2.3.3 of Volume II. Parameter values for the SF— 121 pitch control system are 
given in Section 3.8.3 of Volume II. These parameter values were established 
by the control system analyses described in Section 3.7 of this volume. 

2.8.3 Yaw Control System Model 

The yaw control system model is depicted on figure 2-29. Table 2-3 
gives the parameter values required to implement several common combinations 
of yaw control system types and forward loop controllers. Feedback variable 
selection is controlled by the four selector gains; K. for 0, K r 

r s 8 

for yaw rate, K for \j/ , and Kj r for rj^y; which are functions of CS^ 
and generally assume values of 0. or 1. Those gains which can assume any 
value are indicated in Table 2-3 as being ^ 0 . The exact values of these 
gains must be established by off-line control system analyses. 

Figure 2-29 represents the yaw control system for a VATOL airplane in 
which the cockpit rotates as hover is approached: cockpit control of yew is 

switched from pedals to lateral stick and the yaw trim function switches from 
pedal trim controller to lateral stick trim controller. This feature can be 
eliminated from the model by removing the CS^ W dependence from the cockpit 
controller inputs to the system. 
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Figure 2-28. Generic Pitch Control System *>del 
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+* NA o Not applicable, i.e. the control system type and forward path rontroller are independent of this parameter. 















































The lateral acceleration (n ) to FCS yaw control command 

J'CQ 

(«YAy) feedback was added to decouple a coupled roll-spiral mode which 
appeared in the SF— 121 in the speed range of 60 to 200 kt. (See Section 3.7 
In this volume). The feedback can be removed simply by zeroing the lateral 
acceleration gain (K ). 

y 

Equations for the yaw control system model are given in Section 2.3.4 
of Vblume II. Parameter values for the SF— 121 yaw control system are given in 
Section 3.8.4 of Volume II. These parameter values were established by the 
control system analyses described in Section 3.7 of this volume. 

2.8.4 Heave Control System Model 

The heave control system model is depicted on figure 2-30. It 
receives input rate commands from the heave rate cockpit controller, compares 
these inputs with a heave rate feedback, and adjusts thrust to reduce the 
error according to a proportional (K, » 0.) or proportional plus 

e I 

integral (IU ^0.) forward loop control law. The heave control system 
e I 

operates only at low speed (CS SW =1.). The error gain (K z ) and error 
integral gain (K z ) must be determined by off-line control system 
analyses. 6 1 

Equations for the heave control system model are given in Section 

2.3.5 of Volume II. Parameter values for the SF— 121 heave control system are 
given in Section 3.8.5 of Volume II. These parameter values were established 
by the control system analyses described in Section 3.7 of this volume. 
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Figure 2-30. Generic Heave Control System Model 



2.9 Relative Orientation of Aircraft Body Axes and Inertial Axes 

Most simulations Involving vehicle dynamics require a transformation 
to define the orientation of a set of vehicle-fixed body axes relative to an 
earth-fixed or Inertial axis system. There are two basic formulations of this 
transformation: Euler rate equations and direction cosines. Because of their 

simplicity and almost universal applicability, the Euler rate equations are 
used for most aircraft simulations. Their one drawback is a singularity at o 
• 90 degrees. Even in aircraft simulations, such as air to air combat, where 
rapid large angle maneuvers are involved, the Euler rate equations are used by 
introducing approximate continuous forms of the equations when » is close to 
90 degrees. These approximations work well and have tolerably small errors as 
long as o • 90 degrees occurs only momentarily. It would be impossible to 
trim the aircraft at o - 90 degrees with the approximations present. VATOL 
aircraft routinely maneuver for long periods with o close to 90 degrees. In 
addition they must be trimmed near 90 degrees for initiating hover analyses. 
For these reasons, the direction cosines formula! '"'•n was adopted for the VATOL 
simulation math model since it is continuous everywhere. Because nine (vice 
three Euler rate equations) equations are involved, a slightly increased 
computation load is imposed on the simulation computer. Reference (j) 
quantified this increase to be 0.04 millisecond per integration using an Adams 
2nd order Integration algorithm; the Euler rate equations useo 70 
milliseconds per Integration while the direction cosines used 0.74 
milliseconds. In real time manned simulation activities to date using this 
VATOL model and requiring a frame time of approximately 50 milliseconds, this 
increased load has not been noticed. Neither has it noticeably Increased 
expected rur times in applications of VATLAS to generate dynamic check cases 
for the manned simulation. 

Alternatives to the direction cosines formulation which have not been 
explored are as follows: 

1. Rotated Inertial Axis System - Choose an Inertial axis system in 

which the Xj and Zj axes are rotated relative to the standard 
Inertial system in which the Zj axis is aligned with gravity. 
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The Euler rates and angles would then be applied to the rotated 
axis system. For example, figure 2-31 shows that, with the 
inertial axis system rotated 45 degrees, VATOL hover occurs at 

- 45 degrees (where or is pitch angle referenced to the 
rotated axes) and the singularity will not occur until the 
aircraft rotates another 45 degrees off the vertical. Similarly, 
straight and level flight in the rotated axes will occur with 
®r « -45 degrees. 

2. Qi Uernions - Computational speed and universal applicability are 
cued as advantages of quaternions. A limited investigation by 
Vought of these claims (ref ( j ) ) found that quaternions are 
slightly faster than both Euler rate equations and direction 
cosines (0.66 milliseconds per iteration vs 0.70 milliseconds for 
Euler rate equations and 0.74 milliseconds for direction cosines) 
and are continuous everywhere. This speed improvement is "down 
in the noise" of real time manned simulations which operate in 
the 40 to 60 millisecond range. 

The equations for the direction cosine formulation are detailed along 
with the aircraft equations of motion in Section 2.4 of Volume II. This 
formula* ion incorporates relations developed in reference (k) which maintain 
the orthonormality of the transformation by correcting the direction cosines 
for integration errors. 
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Figure 2-31. Relation of Standard and Rotated Inertial Axis Systems 
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2.10 Pseudo-Pilot Functions 


To effectively apply the VATOL simulation model in an off-line mode 
required the incorporation of pilot-like (or pseudo-pilot) logic and inputs 
into VATLAS. These pseudo-pilot options include: 

1. Apply open-loop inputs to any combination of cockpit controllers 
using data tables. 

2. Equations for pilot-flown transition 

3. Equations for pilot-flown stationkeeping 

The table input functions can be applied in parallel with pilot-flown 
transitions or stationkeeping. Since the pseudo-pilot functions are not 
required for the VATOL manned simulation, they are not detailed in Volume II. 
Instead they will be presented here. 

2.10.1 Open-Loop Cockpit Controller Inputs 

The equations for the open loop cockpit controller inputs are quite 

simple: 


6 LNGSTK ■ f PPj( t ) 

6 LATSTK ■ fypJ*) 

6 PED * f pp 3 ^) 

Z e c “ f PP 4 ( l ) 

6 THR0T * f pp 5 ( t ) 

where the fpp's are user defined tables 
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*LNGSTK 1s fw*e-aft motion of the control stick 


*LATSTK left-right motion of the control stick 
*PED is Pedals motion 

Z e is heave rate commanded by the heave rate controller 
c 

*THR0T * s throttle motion 

Ordinarily the input functions are arranged to be steps, doublets, pulses, and 
other test inputs but, because of their table input format, can be made 
completely random. VATLAS is arranged so that inputs can be applied to any 
and all of the controllers simultaneously. In addition, the inputr can be 
different for each controller since each controller has its own function table. 

2.10.2 Pseudo-Pilot Transition Equations 

The pilot-flown transition equations require tables of trim pitch 
angles and throttle settings as a function of ground speed. Thus an analysis 
of transition trims is required before the pseudo-pilot can "fly" the 
transition. The equations assume that the pitch control system is rate 
command-attitude hold at high speeds (CS SH . o.) and attitude command at low 
speeds (CSj W . 1.) The equations are given below: 


1. 

V e 


? 7 1 

v ♦ w‘ 

2. 

• 

v e 

■ (uu + 

vv ♦ ww)/ 

3. 

< v e 

est' ' 

V e ♦ v e aT 

4. 

4 

- f pp 6 

<V 

5. 

®c 

4 

■ H 

( W 
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6« m ( 9 ^ — 9^ ) / AT 

7 * 6 LNGSTK * V 1 • CS S«) + (•Cj cs sw 

8 * 6 THR0T " f PP ? ( v e> 
where V g is inertial speed 

V e is rate of change of inertial speed 

V a is estimated inertial speed at the next update assuming 
e EST* 

V g is constant 

aT is the time between updates 

# » • 

u, v, w are the rates of change of the body axis components of 
inertial speed 

© r is the trim pitch angle at speed V. 

L 1 e 

© r is the trim pitch angle at speed V 0 
l 2 e EST 

« 

©0 is the pitch rate required over th next aT seconds to reach 
at the next update 

u 2 

q r is ‘he value of the input integrator when CS<- U changes 

from 0. to 1.0 (i.e. when the pitch control system switches 
from rate command-attitude hold to attitude command). 

To initiate a transition $ynR0T retar ^ e ^ below trim to start a 
deceleration, then as the speed decreases to 80 to 100 kt, the throttle is 
advanced to the proper trim setting. These throttle changes are controlled by 
the values in f pp (V g ). When CS$ W changes to 1.0, the heave control 
system is engaged 7 and adjusts the thrust level to maintain the preset heave 
rate (normally = 0 ft/sec.). A sample pilot-flown transition is shown in 
figure 3-32 and discussed in Section 3.8.9. 
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2.10.3 Pseudo-Pilot Stationkeeping Equates 

The development of the stationkeeping control equations assumed that 
the aircraft pitch and yaw control systems are attitude hold and that the 
heave rate command system was operative. It also assumed that the three 
inertial position loops are closed by the pilot who has both inertial rate 
and position available. The stationkeeping control laws are depicted in 
figure 2-32. The pseudo-pilot's goal is to maintain aircraft position at 

« Y e * 0 and Z g at some desired altitude. Since Slng$tk contro ^ s Pitch 
attitude which in turn controls fore-aft motion along the Z body axis (recall 

o«s90 degrees while stationkeeping) and « LATSTK controls yaw attitude which 

in turn controls right-left motion along the Y body axis, the X and Y« 

C c 

rate and position data are resolved into the body axis system. Longitudinal 
and lateral stick commands are then made proportional to Z and Y body axis 
errors respectively. Similarly, the heave rate commands are made proportional 
to heave and heave rate errors or, equivalently (with e * 90 degrees), X body 

axis position and rate. By reference to figure 2-32, the pilot-flown 

stationkeeping equations can be written as follows: 

t 

U \ - - d 13* x e * We> - d 23< Y e * W ~ - \ «"•» - V 

0 

t 

2. Y . _d 12 (x e + K v X e ) -d 22 (Y e + K v Y e ) ? - | vdt - K v v 
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t 



4 * 6 LNGSTK - K Z a Zp 
e r S 
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where 


x e * Y e , Z e are the error quantities along the X, Y, Z 
s bo£y axis assuming o - 90 degrees. The position contributions 
(i.e. }wdt, fvdt, (udt) in the expressions at the far right of 
equations 1, 2, and 3 are exact only if d^* d 23 , d^ 2 , and 
d 22 are constants (i.e. the relative orientation of the body 
and inertial axes does not change). 

K u . K v , are the rate feedback gains and must be determined by 
off-line control system analyses. 

K X » K y , K 7 are the error gains and must also 
e s U e s 

be determined by off-line control system analyses. 

A sample pilot-flown turn over a spot in a 35 kt wind at a commanded 
20 degrees/sec turn rate is shown in figure 3-33 and discussed in Section 
3.8.10. This maneuver combined the pseudo-pilot stationkeeping equations with 
an open loop pedal input. 
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3.0 APPLICATION OF VATLAS TO THE VOUGHT SF-121 AIRPLANE 


This section describes the application of VATLAS to develop and 
demonstrate a FCS for terminal operations of the Vought SF-121 airplane. The 
general FCS design procedure depicted on figure 3-1 was followed for this 
application, ina primary outputs of this iterative process are specification 
of actuator inputs and FCS laws and gains. Required inputs include aircraft 
data, operational conditions, desired FCS configuration, and inputs for 
testing FCS performance. Also required are an initial set of actuator input 
specifications and FCS laws and gains. Each step in the procedure (denoted by 
the blocks in figure 3-1) is discussed in detail below. Cal cions are 

performed by VATLAS, which is used to generate 6D0F nonlinec -.rs and 

response time histories to test inputs, and LINA, Vought 's li analysis 
routine. 

The purpose of this VATLAS application to the SF-121 was to develop a 
baseline FCS which would enable a pilot to perform transition and hover flight 
tasks in a moving base simulation of the airplane. This FCS is to be capable 
of trimming the aircraft anywhere in the transition corridor and of meeting 
the maneuver requirements of MIL-F-83300 and AGARD 577. There ?ce no 
requirements that this FCS be optimized with regard to pilot workload, 
augmentation level required to accomplish terminal operations, ride qualities, 
and similar items which can be considered when defining a baseline FCS but can 
not be completely evaluated without simulator studies and flight tests. 
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Figure 3-1. Flight Control System Design Procedure 







;.l SF-121 Airplane Data 


The data required by VATLAS to model an airplane have been discussed In 
Secton 2.0 of this volume. Requisite SF-121 data have been Incorporated In 
Volume II of this report as follows: 

o Mass Properties Data - Section 3.1 
o Geometry Data - Section 3.2 
o Aerodynamic Data - Section 3.3 
o Propulsion System Data - Section 3.4 
o Inlet Ram Forces and Moments Data - Section 3.5 

o Reaction Control System Data - Section 3.6 

o Actuation System Data - Section 3.7 

The mass properties data in Volume II represent only one SF-121 loading - 
design mission stores (2 Sidewinder and 2 Sparrow missiles and gun plus 400 
rounds of ammunition), gear down, 1000 lb. fuel. This is a landing condition 
and served as the loading for the baseline FCS design. Tables 3-1 through 
3-4, taken from reference (c), provide more extensive mass properties data for 
design mission stores on and off and gear up and gear down. 

SF-121 control system data are given in Section 3.8 of Volume II. These 
data are outputs of the FCS design procedure; they were not required or 
available at initiation pf the design process described herein. 



SF— 121 MASS AND INERTIA DATA: DESIGN MISSION STORES, GEAR UP 
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TABLE 3-2 
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SF-121 MASS AND INERTIA DATA DESIGN MISSION STORES OFF, GEAR DOWN 
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3.2 SF— 121 Operational Conditions 


As stated above, the baseline FCS development was to consider only 
terminal operations of the SF-121. Terminal operations Include approach to 
transition, transition or reconversion, and hover flight conditions. It Is 
assumed that these operations occur at low altitude In a speed range of 0 to 
200 kt. At 200 kt, the SF-121 cruises In aerodynamical ly supported 
(conventional) flight at sS9.S degrees while at hover the airplane is fully 
thrust supported at a&90 deg. The specific operational conditions selected 
for linear analysis and design of the baseline FCS were V A . o, 10, 40, 60, 

80, 120, and 200 kt. In level flight at VL (vertical Landing) weight and 
Inertias. The closely spaced conditions at V A » 40 , 60, and 80 kt emphasize 
that portion of the SF-121 transition where the aircraft trim characteristics 
change most rapidly with speed. (See Section 3.6 for the SF-121 trim data.) 
Landing was considered more critical to the design since VL trim thrust (thus 
control power) Is lower than VTO (Vertical Takeoff) thrust and VL is a more 
demanding piloting task than VTO. 

Another operational condition imposed by the NASA-Ames S.01 simulator 
limitation on cockpit pitch angle but certainly within the realm of VATOL 
possibilities was that. In hover, the cockpit centerline will be approximately 
perpendicular to the fuselage centerline; l.e. the pilot looks out the 
airplane Z body axis during hover. This cockpit (pilot) orientation is 
similar to that of the “nutcracker" aircraft concept. To maintain harmony 
between cockpit controller Inputs and pilot - perceived aircraft motions, the 
baseline FCS has to provide for swapping pedals and lateral stick functions 
during transition. Thus, lateral stick always commands cockpit roll and 
pedals always command cockpit yaw. 
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3.3 Baseline FCS Configuration Selection 


It was immediately obvious that at least two modes were required for the 
SF-121 baseline FCS configuration; one mode for hover or thrust-supported 
flight and one for conventional or aerodynamically-supported flight. A switch 
point or blending region between these modes has to be established. Sections 
3.6 and 3.7 detail the development of this mode changing. 

The next step in FCS configuration selection was to consider how control 
forces and moments are produced for the SF-121 airplane. Table 3-5 summarizes 
this information for each of the airplane six degrees of freedom and for hover 
and conventional flight modes. The table specifies the primary force or 
moment generator and cockpit controller, if available, for each degree of 
freedom. The cockpit heave controller in hover flight is a separate 
controller if closed loor heave or heave rate control is adopted; otherwise it 
is the throttle. With the exception of the role reversal of pedals and 
left-right stick motion and the possible addition of a heave controller in 
hover flight, the functions of the cockpit controllers are fairly standard. 

The dual-mode requirement for the FCS and similarity in number and func- 
tion of the cockpit controllers led to adoption tor the SF-121 of the augmen- 
tation levels and functions of the baseline FCS of the lift/cruise fan V/STOL 
airplane studied in reference (e). The selection of these levels and func- 
tions for the lift/cruise fan airplane FCS was based on an extensive founda- 
tion of analysis, simulation, anu flight testing of terminal operations of 
V/STOL airplanes. These levels and functions were demonstrated in reference 
(e) to be adequate for launch and recovery operations on a DD963 type ship. 
There was no reason to doubt that the same FCS concept would suffice as a 
baseline FCS for terminal operations the SF-121 airplane. 

The SF-121 baseline FCS configuration is summarized by degrees of freedom 
controlled, system type, cockpit controller, and flight mode in Table 3-6. 

The only deviation from the lift/cruise fan V/SfOt baseline FCS is the 
substitution of a rate cornnand-attitude hold system for an attitude command 
system in the pitch axis for the conventional flight mode. The pitch angle 
changes required during VATOL transition are large (50-60 degrees). With an 
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CONTROL FORCE AND MOMENT GENERATION ON THE SF-1Z1 
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SF-121 BASELINE FCS CONFIGURATION 
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attitude command system limited to 15 or 20 degrees pitch command authority 
for the stick, the pilot has to spend considerable time on the trim button to 
keep the stick in a desirable range of travel. The rate command-attitude hold 
system eliminates the requirement for trim button Inputs and maintains the 
attitude stability supplied by the attitude command system. 
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3.4 Test Cases to Evaluate PCS Performance 


The test cases selected to demonstrate and evaluate FCS performance are as 
follows: 

o Steps, doublets, and/or pulses imposed on various cockpit controllers 
in both hover and conventional modes and in the mode switching f 1 J jht 
regime. 

o Pseudo-pilot flown transition 

o Pseudo-pilot flown turn over a spot in a 35 kt wind 

Since these runs will be produced by the full nonlinear simulation program 
(VATLAS), they should highlight any problems due to actuator position or rate 
limits, nonlinear aerodynamics, thrust dynamics, kinematic coupling, and 
similar nonlinearities. The pseudo-pilot transition and stationkeeping runs 
will indicate potential piloting problems in attempting to perform transition 
and hover tasks. 
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3.5 Actuator Input Specification 


The SF— 121 has fourteen actuators which are to be linked with four FCS 
variables. The FCS variables are commanded thrust and roll, pitch, and yaw 
control commands. The roll, pitch, and yaw commands are normalized such that 
♦1. commands all available control power In the appropriate direction. The 
fourteen actuators are as follows: 

o Left and right wing trailing edge flaps (elevons) 
o Vertical stabilizer control surface (rudder) 

o Left and right engines pitch thrust deflection 

o Left and right engines yaw thrust deflection 

o Left and right engines thrust (l.e. propulsion system) 

o Left and right wing tip RCS jets 

o Canard trailing edge flap 

o Left and right wing leading edge flaps 

The inputs to the left and right wing leading edge flap actuators and to 
the trailing edge flap of the canard are specified in reference (c) to be 
functions of angle of attack. Inputs for each of the other actuators except 
the propulsion system were formed by multiplying the appropriate normalized 
roll, pitch, or yaw command by the maximum deflection controlled by that 
actuator. For example, the normalized pitch control ( 6 p xtch^ was ^ nked *° 
the right and left wing elevon actuators (which are the primary generator of 

pitch control power during conventional flight) and to the left and right 

engines pitch thrust deflection actutators (which are the primary generator of 
pitch control power during hover flight). All actuators were assumed to be 
driven full time. The gains for the elevon actuator inputs were 25 degrees, 
the maximum symmetric elevon deflection, and the gains for the pitch thrust 
deflection actuators were 15 degrees, the maximum pitch thrust deflection. 

The actuator inputs are depicted in figure 3-2. This arrangement is the 
initial and, as will be shown, the only actuator input specification required. 
Note that, In addition to the basic links, the arrangement anticipates use of 
differential left-right pitch thrust deflection for roll control (gain 
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« r ) and a crossfeed of roll control to rudder (gain K 4 « r ) and 
yai thrust deflection (gain « r ). As will be shown, th^se additions 


were not required. 
Volume II. 


The actuator input gains are specified in Section 3.7 


.1 of 
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3.6 Trim Characteristics Evaluation 


Once the actuator Inputs were specified, the airplane could be trimmed 
using FCS control variables. Two types of trim were calculated by VATIAS: 
hover trims in a 35 kt wind whose direction varied from head wind to tail wind 
and transition trias. Transition trias assuaed straight (# « 0) and level 
(y * 0) flight and were calculated for three deceleration levels (V A * o, 

-9.1, -0.2 g) along the flight path for a speed range of 0 to 200 kt. The 
primary transition trim variables - pitch angle, thrust level, and normal ed 
pitch control ( itch) “ are shown on figure 3-3. Similarly the primary 
hover trim variables - pitch and bank angles, thrust level, angles of attack 
and sideslip, and normalized roll, pitch, and yaw controls - are shown on 
figure 3-4. 

Figures 3-3 and 3-4 indicate that the SF-121 has sufficient control power 
for trim in 35 kt winds and along reasonable reference trajectories in the 
airplane's transition corridor. The maximum percentage use of available pitch 
control power to trim is 74%which occurs at 80 kt in a 0.2g decelerating 
transition. The maximum percentage use of available roll control power to 
trim is 9%which occurs in a 35 kt crosswind - -90 degrees) while 

that for yaw control power is 17% which occurs in a 35 kt wind oriented 75 deg 
to port or starboard. 

Figure 3-5 compares pitch, roll, and yaw control power available with trim 
(for a 0.2g decelerating transition) and flying qualities specification 
requirements. Because o* the rotated SF-121 cockpit in hover the flying 
qualities roll (yaw) control power requirements are compared with SF-121 yaw 
(roll) control power. AGARD 577 requirements are assumed to be applicable 
from C to 35 kt, which is the range for hover and low speed flight specified 
in MIL-F-83300. The AGARD 577 requirements indicated on figure 3-5 show the 
control power range designated to be typical for attitude stabilized V/STOL 
aircraft frr maneuvering, trim, and gust regulation functions. The control 
power available for maneuvering and gust regulation is also compared with the 
flying qualities requirements on figure 3-5. As can be seen the SF-121 has 
adequate control power for meeting the trim and maneuvering specification 
requirements with some margin left for gust regulation. 
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Figure 3-3. SF— 121 Trim Requirements for Level Flight Decelerations (Sheet 2 of 2) 







Figure 3-4. SF-121 Trim Requirements for Stationkeeping In a 35 kt Wind (Sheet 1 of 4) 
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Figure 3-4. SF-121 Trim Requirements for Stationkeeping in a 35 kt Wind (Sheet 2 of 4) 
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Figure 3-4. SF— 121 Trim Requirements for Stationkeeping in a 35 kt Wind (Sheet 4 of 4) 
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The results of the trim characteristics evaluation have shown that the 
SF— 121 has sufficient control power and that the actuator Input specification 
of Section 3.5 is adequate. Thus no iterations were required of the actuator 
input specification design loop shown in figure 3-1. A more comprehensive 
trim analysis might have indicated a nonlinear blending of the elevon and 
pitch thrust deflections is better than the linear blending which has been 
adopted. For example, the elevon could be programmed to provide all the trim 
control power to as low speeds as possible. This would save thrust deflection 
for the maneuvering and gust regulation functions. 

Another result of the trim characteristics evaluation was the selection of 
60 kt as the switch point between the hover and conventional modes of the 
FCS. To obtain this result the transition trim data of figures 3-3 and 3-4 
were examined using the following criteria: 

o More than 50%of the aircraft weight supported by thrust (i.e. 

transition more than 50^complete) . This occurs at approximately 65 
kt. 

o Since surge and sway control in hover are generically the same for 
the SF-121, similar capabilities and characteristics are desirable in 
the paired pitch/surge and yaw/sway degrees of freedom. Thus it was 
decided that pitch control margin and yaw control power should be 
approximately equal at and below the switch point. As shown by 
figure 3-5 this criterion is first met near 60 kt when the two 
control powers equal 0.7 rad/sec^. 

The final result of the trim characteristics evaluation was the provision 
of dimensional stability derivatives for developing the FCS laws and gains. 
These derivatives are an important by-product of VATLAS trim calculations; 
they are automatically generated at each trim point. Tables 3-7 and 3-8 
summarize SF-121 derivatives at the speeds selected for FCS development. The 
derivatives are based on a transition trim along an unaccelerated straight and 
level flight path. 
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SF-121 LONGITUDINAL STABILITY DERIVATIVES AT FLIGHT 
CONTROL SYSTEM DESIGN CONDITIONS 
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SF-121 LATERAL DIRECTIONAL STABILITY DERIVATIVES AT FLIGHT 
CONTROL SYSTEM OESIGN CONDITIONS 
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3.7 Specification of Baseline FCS Laws and Gains 


Linear analysis, supported by LINA, was applied to quantify the baseline 

FCS conflguraton described In Section 3.3 and Table 3-6. Design analyses were 

performed at each of the speeds specified In Tables 3-7 and 3-8. The first 

step In these analyses was to calculate the unaugmented or bare airframe 

transfer functions using the stability derivatives In the tables. These 

airframe functions are given In Tables 3-9, 3-10, and 3-11. 

« 

Longitudinal (pitch and heave) FCS development will be discussed In 
Section 3.7.1 followed by lateral directional (roll and yaw) development In 
Section 3.7.2. Discussion of the hover FCS design details are focused on the 
V A > 10 kt case while conventional FCS details are focused on the V A - 120 
kt case. The FCS for other speeds differs only in gain from the 
representative 10 and 120 kt cases. 

Before the baseline FCS configuration (Table 3-6) could be quantified, 
several issues having general applicability to the design process had to be 
resolved. These were: 

o How should FCS gains be scheduled 

o What attitude signal should be fedback for attitude stabilization in 
the hover mode 

o What are desirable bandwldths for the attitude stabilization and 
heave rate loop.. 

The first issue was resolved In favor of gain scheduling as u function of 
airspeed. Pitch angle scheduling was considered but quickly eliminated; a 
wide range of pitch attitudes can be attained at any speed, thus a high 
probability of Incompatible FCS gains exists with possibly disastrous 
consequences. 

The Issue of attitude feedback signals for hover stems from the discon- 
tinuity at o ■ 90 degrees In the Euler angle transformation. If the aircraft 
body axes are aligned with standard North-oriented Inertial axes at 
o . p • y m 0 deg. then 4 and yt'must change by 180 degrees to maintain the 
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34.37 (- 6 . 00 ) (-. 01 28 ) ( 5 . 32 ) 11.68 (- 4 . 329 ) ( 3 . 631 )(. 325 ) 4.26 ( 2 . 264 ) (- 3 . 1 51 )(. 150 ) 
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First number in box is root locus gain of tranfser function 




















SF— 121 UNAUGMENTEO AIRCRAFT LONGITUDINAL TRANSFER FUNCTIONS 
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First number In box is root locus gain of tranfser function 






















SF-I21 UNAUGMENTED AIRCRAFT LATERAL DIRECTIONAL TRANSFER FUNCTIONS (V A . 60, 80, 120 and 200 kt) 



means 1st order root at 




































































SF-1ZI UNAUGMENTED AIRCRAFT LATERAL 01 RECTI ON AL TRANSFER FUNCTIONS (V A . 60, 40. 10 and 0 kt) 
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First number In box Is root locus gain of transfer function 






























consistency of the transformation as the aircraft Xg axis goes through 
vertical. These large angle changes would be imposed in the appropriate 
feedback paths and produce cockpit motions which are disorienting to the 
pilot: For example, to slow down in surge the pilot simply pulls back on the 

stick to bring the aircraft nose through vertical. He doesn't command or 
expect cockpit roll or yaw motions during this maneuver. With Euler angle 
feedbacks based on the standard North-oriented inertial axes the pilot will 
experience uncommanded cockpit yaw and roll motions which (if control isn't 
lost and/or the pilot doesn't punch out) will place him in the orientation he 
expects. Analysis of this problem indicated that uncommanded cockpit motions 
could be avoided in hover by feeding back the integrals of body axis roll, 
pitch, and yaw rates as attitude feedbacks. (Also, the rotated inertial axis 
system discussed in Section 2.9 could be applied.) Thus the attitude feedback 
issue was resolved as follows: feedback the integrals of body rates in the 

FCS hover mode and standard Euler angles in the FCS conventional mode. In 
addition, the body rate integrators must be initialized to equal Euler angles 
upon switching to hover mode. Similarly, upon switching to conventional mode, 
the trim systems must be reinitialized to avoid discontinuities in attitude 
feedbacks. 

Desired attitude loop bandwidth was set at 3.0 rad/sec. This selection 
was based on experimental evidence (references (1) and (m)) that pilots strive 
for closed loop bandwidths of approximately 3.0 rad/sec when they manually 
close the loop. This provides adequate separation between the flight path and 
attitude modes of the aircraft and allows the pilot to attain reasonable 
manual closures of the flight path control loops. In addition to reducing 
pilot workload, a 3.0 rad/sec attitude bandwidth also provides good regulation 
against distrubances. 

Desired heave rate bandwidth was set at 1.5 rad/sec. Similar to the 
attitude bandwidth selection, this rate bandwidth allows the pilot to easily 
close the manual heave position loop in the desired bandwidth range of 0.7 to 
1 .0 rad /sec. 
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An issue of lesser importance to quantifying the baseline FCS was the 
specification of parameter values for the pilot input shaping networks in the 
rate command-attitude hold and attitude command systems. Even though these 
parameter values will be established via the piloted simulator studies, 
initial "guess" values are required at this point in the design process. 
Therefore the first order lag in the attitude command system input shaping 
network was set at 4.0 rad/sec in accordance with Vought experience while the 
breakpoint of the lead term in the rate command-attitude hold systems input 
shaping network was set at 3.0 rad/sec which was used in reference (e). Also 
in accordance with reference (e) input command authority was set to +15 deg. 
for the attitude command systems and +20 deg/sec for the rate command-attitude 
hold systems. 

3.7.1 Longitudinal FCS Development 

The baseline configuration for the pitch control system is shown in figure 
3-6. The values adopted for many of the generic pitch control system 
parameters are evident in a comparison of figures 2-28 and 3-6. The baseline 
heave control system is essentially the same as the generic system shown in 
figure 2-30. "tost features of these configurations have already been 
discussed in Sections 3.3 and 3.7. The functional descriptions which follow 
are intended as a review but provide additional details where required. 

The forward path of the pitch FCS consists of a straight gain (K q ) and 
a parallel integrator (gain ■ K q ) that provides automatic trim and 6 

e Il 

boosts the low frequency gain of the loop. K q is programmed as a func- 
tion of airspeed. The feedbacks in the hover fiode include pitch rate through 

a gain (K q ) i n parallel with the integral of pitch rate (Qint)> * n the 

conventional mode the feedbacks include pitch rate through K in parallel 

q q 

with pitch angle (®). The pitch FCS has two pilot inputs, the longitudinal 

stick trim button (« q ) an d longitudinal stick deflection UlnGSTk)* 

In the hover mode, ’ s in P ut through a gain (K q ) and first order 

shaping network (time constant *^puc H ). In the conventional mode, 

5 LNGSTK i n P ut - through K g and a parallel integrator (gain • K fl . 

x H c./» 
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Figure 3-6. ST-121 Baseline Pitch Flight Control System 








TpiTCH is zeroed. Thus, \NGSTK COfl,m * nd s are proportional to the integral 
of pitch rate (qj NT ) in the hover mode and to pitch rate (q) in the 
conventional mode. C Logic for initializing the pitch rate integrator, manual 
trim integrator, and pitch input integrator at the FCS mode switch point is 
based on that indicated on figure 2-28. 

The heave FCS is operative only in the hover mode. In conventional mode 
the pilot controls thrust directly with manual throttle. In hover mode he 
normally controls heave (and thus thrust) through his heave rate command 
lever. Manual throttle is available during hover but its use would likely be 
limited to emergency conditions requiring sudden large changes in thrust. 
Similar to the pitch FCS, the forward path of the heave FCS has a straight 

gain (K^ ) j n parallel with an integrator (gain = K z ) for automatic 
e e i 

trim ard low frequency gain boost. The feedback is heave rate which is 
compared with the heave rate command lever inputs to form the forward path 
error signal. It will be demonstrated in Section 3.8 that the performance of 
both modes of the pitch and heave FCS laws and gains developed here is 
adequate. Thus there was no need to iterate the laws and gains loop of the 
FCS design procedure (figure 3-1). Complete specification of the SF-121 
baseline pitch and heave FCS gains is given in Sections 3.8.3 and 3.8.5, 
respectively, of Volume II. 

3.7. 1.1 Hover Mode System Analysis 

A root locus sketch of the pitch FCS hover mode loop closure at = 10 
kt as a function of K is shown in figure 3-7. The zero at s = -0.2 = 

- K 0 and pole at the origin result from the forward path compensation. 
e Il 

The zero at s * -2.0 * -1/K results from the pitch rate feedback gain 

M q 

and was placed to effect a good closed loop damping ratio (0.7 - 0.8) at the 
desired system bandwidth (approximately 3.0 rad/sec). With the exception of 
the elevon and pitch thrust deflection actuators both having poles at s = -20, 
the other open loop poles and zeros are those of the unaugmented airplane 

q INT/ 6 P ITCH trans ^ er function. As can be seen, the selected K q 
(= 10/rad) provides dominant closed loop roots which appear to provide a 
larger than desired system bandwidth. As indicated by the closed loop 
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Figure 3-7. Root Locus of Pitch FCS Hover Mode Loop Closure 



q INT^INT frequency response (figure 3-9), this is not the case. The 
equivalent closed loop bandwidth is 2.8 rad/sec which is slightly less than 
desired. 

A root locus sketch of the heave FCS loop closure at V A - 10 kt with the 

pitch FCS hover node loop closed is shown in figure 3-8. The forward path 

gain (K z ) varies along the locus. The zero at s • -0.2 « -K z and 
e ej 

pole at the origin result from the forward path compensation. The pole at s * 
-5.0 represents the thrust dynamics at the trim thrust setting. The remaining 
"open loop" pole at s - -.00905 is the only one of the z e / « HEA vg| q INT . q * « PITCH 
transfer function poles not effectively cancelled by zeros. The designation, 

• I • 

Z e^*HEAVE|qiNT» q * *PITCH* frdi cates Z e^*HEAVE transfer function 
which results when the pitch FCS loop is closed (i.e. q INT and q are fedback 

to 6pjj£ H ). The value of K z selected (■ 500 lb/ft/sec) provides a 
heave rate loop bandwidth of e 1.41 rad/sec which is acceptably close to the 

desired 1.5 rad/sec bandwidth. This is demonstrated by the closed loop 

• • 

l Jit, frequency response in figure 3-10. 

C c 

3.7.1 .2 Conventional Node System Analysis 

A root locus sketch of the pitch FCS conventional mode loop closure at 
a 120 kt as a function of is shown in figure 3-11. The zero at 
s - -0.2 « -K and pole at tHe origin result from the forward path 
e Il 

compensation. The zero at s » -2.0 - -1 /K results from the pitch rate 

M q 

feedback gain and was placed to obtain a closed loop bandwidth of 
approximately 3.0 rad/sec. With the exception of the elevon and pitch thrust 
deflection actuators both having poles at s * -20, the other open loop poles 
and zeros are those of the unaugmented airplane ®/« pitch tran sfer function. 

Note the unstable short period pole which arises from the static instability 
designed into the SF-121. As confirmed by the ®/® c frequency response 
(figure 3-12), the selected K (» 7. /rad) provides a 3.15 rad/sec attitude 
bandwidth. The e/® c frequency response does not include the effects of the 
longitudinal stick shaping network [(s + 3.33)/s] which makes the conventional 
mode pitch FCS a rate command-attitude hold system. 
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Figure 3-8. Root Locus of Heave FCS Loop Closure With Pitch FCS Hpver Mode Loop Closed 
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Figure 3-9. Frequency Response of Pitch FCS in Hover Node 





















Figure 3-10. Frequency Response of Heave FCS in Hover Mode 
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Figure 3-12. Frequency Response of Pitch FCS in Conventional Mode 






3.7.2 Lateral FCS Development 


Baseline configurations for the roll and yaw control systems are shrwn in 
figures 3-13 and 3-14 respectively. The values adopted for many of the 
generic roll and yaw control system parameters are evident in comparisons of 
figure 2-27 and 3-13 and of figures 2-29 and 3-14. Many features of these 
configurations have been discussed in Sections 3.3 and 3.7. The functional 
descriptions which follow are intended as a review but provide additional 
details where required. 


The forward path of the roll FCS consists of a straight gain (Kp ) and 
a parallel integrator (gain * K ) that provides automatic trim ami 

Pe Il 

boosts the low frequency gain of the loop. K i S programmed as a 
function of airspeed while Kp i< \ function of control system mode. 

e I 


The feedbacks in the hover mode include body axis roil rate through a gain 
<Kp ) in parallel with the integral of body axis roll rate ( P j nj ) • 1° 


the conventional mode the roll FCS feedbacks are stability axis roll rate 
through K in parallel with roll angle (d). 


The forward path of the yaw FCS is identical in function to that of the 

roll FCS except the gains are K an d K r . Similar to the roll 

e e n 

FCS, K is programmed as a function of airspeed while K r is a 
e e Il 

function of control system mode. The feedbacks in the hover mode include body 
axis yaw rate through a gain (K ) in parallel with the integrai of body 
axis yaw rate fr^y). The feedbacks in the conventional mode reflect re- 
quirements imposed by the bare airframe lateral characteristics of the SF— 121 
(Table 3-10). The airplane has a coupled roll-spiral mode and an inadequate 
relation between the numerator of the r s / 4yaw transfer function and the 
Dutch roll mode. The coupled roll-spiral mode is decoupled by a lateral ac- 
cele. ation (a y ) to 6 yaw feedback through a gain (K a ). The latter 
characteristic makes pure r $ feedback ineffective it a yaw damper and re- 
sulted in the adoption of a pseudo - 8 (t - (g /V A ) $ - r s ) to 6 yaw feed- 
back. The use of s as a yaw damper as well as a turn coordinator is described 
in reference (e). 
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The lateral PCS has four pilot controller inputs - lateral stick trim 

button (a ) pedal trim (« r ), lateral stick deflection 

HBut 'BUT 

( ®LATSTK ) » and deflection (« PE q). The roles of these inputs are a 

function of FCS mode: In the hover mode, the roll FCS receives inputs from 

*P £0 and while the yaw FCS receives inputs from 4 lajstk and 

4 n • ^ the conventional mode, the controller roles are reversed - 

BUT 

®PED and a r become yaw FCS inputs and «latSTK and 6 p n become 

roll FCS in^ifts. The roll FCS is rate command-attitude^Wild in both modes. 

Thus, «p£Q in hover and «i A TSTK ’ n conventional are interfaced with the 
roll FCS through a gain (Kp ) an d a parallel integrator (gain = 

k d ). The yaw FCS in conventional mode is a pseudo-B command system. 

C I 

The pedals are thus interfaced with the yaw FCS through a gain (K r ) a nd 

input filter (time constant Note that, except for deliberate 

sideslips, the pilot should not have to use pedals in the coventional mode. 

In hover mode the yaw FCS is a r INJ command system; 6 LATSTK is 

interfaced through K r an d the input filter. 

c 

Logic for initializing the yaw and roll rate integrators and the trim 
inputs at the FCS mode switch point is based on that indicated on figures 2-27 
and 2-29. This logic helps to smooth the mode change transients. Its 
effectiveness is demonstrated in Section 3.8. 

It will also be shown in Section 3.8 that the performance of both modes of 
the yaw and roll FCS laws and gains developed here is adequate. Thus there 
was no need to iterate the laws and gains loop of the FCS design procedure 
(figure 3-1). Complete specification of the SF— 121 baseline yaw and roll FCS 
gains is given in Sections 3.8.2 and 3.8.4, respectively, of Volume II. 

3. 7. 2.1 Hover Mode System Analysis 

A root locus sketch of the yaw FCS hover mode loop closure at V ft = io kt 
is shown in figure 3-15. The forward path gain (K r ) varies along the 
locus. The zero at s - -0.8 * -K r and one of the poles at the origin 

e Il 
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Figure 3-15. Root Locus of Yaw FCS Hover Mode Loop Closure 



arise from the forward path compensation. The zero at s « -2.0 « -1/K f 
results from the body axis yaw rate feedback and was placed to provide jjood 
loop closure properties. The actuator pole at s * -20. represents the 
dynamics of the rudder and yaw thrust deflection actuators. The remaining 
open loop poles and zeros are those of the unagumented airplane r j^j/ay/^ 
transfer function. The selected closed loop gain (K r , 13 /rad) in 
combination with the roll FCS loop closure (depicted fn figure 3-16) is 
demonstrated by the r* j / r j ^ j frequency response (figure 3-17) to 
produce an attitude bandwidth C of 3.24 rad/sec. 

The root locus of the roll FCS hover mode closure at = 10 kt with the 
yaw FCS loop closed is shown in figure 3-16. Equivalently to the yaw loop, 
the roll loop forward path gain (K^ ) varies along the locus; the zero at 
s = -0.8 « -K and the pole at tfie origin arise from the forward path 
e Il 

compensation; and the zero at s = -2.0 = -1/K results from the body axis 

” P 

roll rate feedback. The remaining "open loop' poles and zeros are those of 

the P int /6 RQ(_I_J transfer function which do not 

r I NT * r * 6 YAW 

effectively cancel each other. Actuator dynamics are not indicated in the 
locus sketch due to their negligible influence. The selected closed loop gain 

(K p = 6/rad) is demonstrated by the P j NT^Pl NT frequency response 
(figure 3-18) to produce an attitude bandwidth c of 3.32 rad/sec. This 
frequency response plot does not include the effects of the pedals shaping 
network [(s + 3 .33 ) / s] which makes the hover mode roll FCS a rate 
command-attitude hold system. 

3. 7. 2. 2 Conventional Mode System Analysis 

As mentioned above, the SF-121 has a coupled roll-spiral mode at the 
airspeeds considered in this FCS development. Traditional design guidelines 
for lateral control law development (also the flying qualities specification) 
prohibit this coupling. The first task in the lateral FCS conventional mode 
development therefore was to decouple the roll and spiral modes. This was 
ar*ompl ished by a lateral acceleration (a ) to « Y AW feedback. Figure 3-19 
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Figure 3-16. Root Locus of Roll FCS Hover Mode Loop Closure With Yaw FCS Loop Closed 



Figure 3-17. Frequency Response of Yaw FCS in Hover Mode 











Figure 3-18. Frequency Response of Roll FCS in Hover 
























figure 3-19. Root Locus of Lateral Acceleration to 6y^ w Loop Closure (V A * 120 kt) 
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shows a root locus sketch of this loop closure at V A . 120 kt. The lateral 

acceleration gain (K a ) varies along the locus. Note the formation of 
the decoupled equivalent roll and spiral modes and the characteristic decrease 
in Dutch roll damping which accompanies lateral acceleration feedback. 

Also mentioned above was the application of psuedo-l ($ « (g/v A ) * - r s ) 
feedback for increasing Dutch roll damping and providing turn coordination. 
Traditionally these functions are provided by a washed out stability axis yaw 
rate feedback. The effectiveness of this feedback depends on the relative 
location of the complex zeros of the r s /« transfer function numerator 
and the Dutch roll mode poles. If these zeros and poles are near the ju axis 
or in the right half plane or the frequency of the zeros is greater than 0.4 
to 0.5 that of the poles, the yaw rate feedback will be ineffective. Both 
these detracting conditions are displayed by the SF-121. Reference (e) 
develops and demonstrates how pseudo-B feedback is an appropriate substitute 
for washed out yaw rate. 


Figure 3-20 provides a root locus sketch of the b to « YAW loop closure 
at 120 kt with the a y to « YAW loop closed. The forward path gain (K r ) 
varies along the locus. The zero at s - -0.2 * -K an d the pole at 

e Il 

the origin are contributed by the forward path compensation. The remaining 
open loop zeros are the zeros of the t/« YAH transfer function while the 
remaining poles are the equivalent lateral directional modes created by the 

to « YA y loop closure. Note that the Dutch roll damping is increased 
and that an unstable low frequency mode is developed by this closure. The . 
unstable mode will be stabilized by the roll FCS loop closure. As is 
demonstrated in figure 3-22, the gain selected for s' loop closure (K r = 

1.1 /rad/sec) in combination with the roll FCS loop closure (sketched fn figure 
3-21) produces a body axis lateral velocity (v 1 fc V A a) to pedal (or t c ) 
frequency response which is similar to that of an unaugmented airplane having 
uncoupled roll and spiral modes and a well damped Dutch roll made. 

The final lateral conventional mode loop closure, the roll FCS, is 
depicted by the root locus sketch of figure 3-21 for V A * 120 kt. The 
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Figure 3-20. Root Locus of Pseudo-e to « YAW Loop Closure With a y to « YAW Loop 
Closed (V A = 120 kt) 
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Figure 3-21. Root Locus of Roll FCS Conventional Mode Loop Closure With a to « v . u and £ to « v . u Loops Closed 














forward path gain (Kp ) varies along the locus. The zero at s • -0.2 * 
-Kp and the pole at the origin are contributed by the forward path 
e U 


compensation. The zero at s 


-3.0 



results from the stability 


axis roll rate feedback »no was placed to effect a reasonable loop closure. 
The actuator pole at s « -20 represents the elevon actuators dynamics. The 
effects of the dynamics of the RCS and its actuators are negligible to this 
loop closure. The remaining "open loop" poles and zeros are those of the 

^/®R0LL* a » ^ * 6 yaw trans ^ er function. As is demonstrated by the 
d/d c frequency reponse (figure 3-23), the gain selected for roll FCS loop 
closure (K p 3 4/rad) produces an attitude bandwidth of 2.7 . rad/sec which 
is slightly e less than the desired 3.0 rad/sec. This frequency response plot 
does not include the effects of the lateral stick input shaping network [(s + 
3.33)/s] which makes the conventional mode roll FCS a rate command-attitude 
hold system. 
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Figure 3-23. Frequency Response of Roll FCS in Conventional Mode 










3.8 FCS Performance Evaluation 


The control laws and gains established by the analyses described in 
Section 3.7 were incorporated into VATLAS. Nonlinear time reponses to various 
cockpit controller inputs and pseudo-pilot flown scenarios were then 
calculated to evaluate the FCS performance. Ten cases selected for evaluation 
and demonstration of this performance are as follows: 

1. Ll jitudinal Stic* doublet at V A - 120 kt commanding 10 deg/sec 
pitch rate for 0.7 sec followed by -10 deg/sec for 0.7 sec (figure 
3-24) 

2. Lateral stick pulse at V A . 120 kt commanding 20 deg/sec roll rate 
for 1.4 sec (figure 3-25) 

3. Pedal step at V A » 120 kt commanding 4 deg/sec of pseudo b (figure 
3-26) 

4. Longitudinal stick doublet at V A ^ 10 kt commanding 10 deg of 

^INT * or 0*7 sec followed by -10 deg for 0.7 sec (figure 3-27) 

5. Lateral stick doublet at V A E 10 kt commanding 10 deg of r j for 
0.7 sec followed by -10 deg for 0.7 sec (figure 3-28) 

6. Pedal pulse at V A = 10 kt commanding 20 deg/sec roll rate for 1.4 
sec (figure 3-29) 

7. Heav rate controller doublet at V A , 10 kt commanding -10 ft/sec 
heave rate for 3.0 sec followed by 10 ft/sec for 3.0 sec (figure 3-30) 

8. Mode switching transient. Aircraft trimmed in a 0.1 g decelerating 
0.5 mile turn at V A = 62 kt, then simultaneous longitudinal and 
lateral stick pulses of 0.7 sec duration commanding 10 deg/sec pitch 
and roll rates are input 0.3 sec into the run (figure 3-31) 

9. Psuedo-pilot flown transition initiated at V A = 200 kt (figure 3-32) 

10. Pseudo-pilot flown turn over a spot (i.e. stationkeeping) m a 35 kt 

wind (figure 3-33) 

Each case is discussed in some detail in Sections 3.8.1 to 3.8.10. 

The doublet, step, and pulse inputs for cases 1 to 7 were givtn sufficient 
mag itudo to induce control system saturation and/or exercise other signifi- 
cant system nonlinearities. As such they simulate fairly vigorous pilot use 


146 



of the SF-121 and its control system. Other than the fact that 0.7 sec is 
approximately twice the time constant of the attitude loops there is no 
rationale for selecting multiples of 0.7 sec for pulse and doublet lengths. 
Similarly 3.0 sec for the heave rate doublet input is approximately tour times 
the time constant of the heave rate loop. Case 8 demonstrates the system 
capability to maintain control in the presence of large simultaneous 
multi-axis pilot inputs at the mode switch speed. Cases 9 and 10 demonstrate 
the capabilities of the VATLAS pseudo-pilot logic and the control system for 
two operational scena-os for the aircraft. 

3.8.1 Case 1 - Longitudinal Stick Doublet at V A , 120 kt 

The pitch rate response is essentially linear until the doublet reverses 
at 1.0 sec (figure 3-24). The linear response overshoots its commanded value 
of 10 deg/sec by approximately 20£ This is consistent with the pitch loop 
frequency response (figure 3-12) which indicates that the system is less than 
critically damped (i.e. it has a 2.5 db peak which indicates a damping ratio 
of approximately 0.4). At 1.0 sec the symmetrical elevon command becomes 
saturated due to the normalized pitch control becoming < -1.0. From this 
point on the response is nonlinear. The resulting response is stable and 
completes the doublet maneuver with considerable lag. Coupling into the 
lateral degrees of freedom is minimal. 

3.8.2 Case 2 - Lateral Stick Pulse at V A = 120 kt 

The lateral stick puls: commands a fairly rapid entry to a 28 deg banked 
turn. A steaoy state turn rate of approximately 4 deg/sec is generated by 
this maneuver. The time history traces shown on figure 3-25 indicate that the 
28 degree ru II angle is attained with no overshoot but reduces a degree or so 
when the lateral stick command is released due to the overshoot in the roll 
rate response. Approximate^ 4 deg of b (adverse yaw) and — 0 . 1 5g lateral 
acceleration are generated during the turn entry. These are rapidly driven to 
near zero (thus coordinating the turn' by the yaw control system when the 
desired roll angle is attained and held. 
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Figure 3-24. Longitudinal Stick Doublet Response at V A - 120 kt (Sheet 3 of 3) 
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Coupling into the longitudinal degrees of freedom is apparent and was 
anticipated. The steady state pitch rate is required by the kinematics of the 
steady turn. The reduced (from trim) steady state pitch angle is due to the 
pltc* control system error returning to zero at steady state; since both pitch 
rate and pitch angle are fed back and steady state pitch rate is non zero with 
a positive value, steady state pitch angle must reduce to keep the pitch 
control system error zero. Increased pitch angle during turn entry is caused 
by the loss of elevon effectiveness at large deflections; turn entry requires 
increased left elevon deflection and reduced right elevon deflection thus 
producing a nose up moment. A nose down moment is produced when the 
differential elevon deflection is removed upon attaining the desired roll 
angle. 

Overall, there are no apparent stability or control problems or unantici- 
pated motions in the aircraft response to the lateral stick pulse. 

3.8.3 Case 3 - Pedal Step at V A = 120 kt 

The pedal step commands the aircraft to perform a wings level (d = 0) 
skidding (a ^ 0) turn (r 4 0). The pedal step on figure 3-?6 is labelled as a 
pseudo-a command input. This is consistent with $ feedback of the yaw control 
system. Since the attitude hold feature of the roll control system keeps 
wings level, the yaw control system feedback becomes yaw rate and the pedal 
becomes a yaw rate controller. 

The a response shows the effects of the two "shelf" v/^ frequency 
response (figure 3-22). The magnitude of the frequency response has a steady 
state level ("shelf") of 57 db which extends to approximately 0.02 rad/sec, 
then the response drops to another "shelf" of 42 db which extends to 
approximately 1.0 rad/sec where the response starts its characteristic drop off 
to -«odb. Thus the a time response should and does demonstrate two distinct 
modes; a shot term response wherein a reaches approximately -4 deg in 2.5 sec 
superposed on a long term response which appears as a drift. The yaw rate 
response has a similar bimodal response. 
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Another interesting aspect of figure 3-26 is the effect of decreased 

elevon effectiveness at large deflections which was discussed in Section 

.8.2. As a increases, the roll moment required to maintain wings level 

increases because of the aircraft dihedral effect (Ci ). This is seen in 

1 $ 

the « R ql|_ trace. The pitch up due to the elevon effectiveness loss must be 
countered by increased nose down pitch moment as shown by the «pjyc H trace. 

Note that the aircraft is nearly out of nose dcvn pitch control («pj TrH 5s _i.o) 
and running out of negative roll control at the end of the 15 sec time history. 
A short time more and the aircraft would be in control saturated situation. 

Thus figure 3-26 demonstrates a e limitation in the transition region of 
flight; a not uncommon characteristic of V/STOL aircraft. 

3.8.4 Case 4 - Longitudinal Stick Doublet at = lo kt 

The longitudinal stick doublet commands the aircraft to pitch through 
vertical tnen H own through trim and then return to trim. Since a. is 
controlled by Ditch attitude the intent of this maneuver is to 

decrease X g then increase it beyond trim X e a.'d then return it to near 
trim. Figure 3-27 shows that the desired pitch ma auver is performed with 
some saturation of the pitch controls; | *p jjCH i exceeds 1*0 and j©jJ is 
limited to 15 deg. The desired effect on X g is not obtained and herein lies 
the most distinguishing flying qualities characteristic of hover control of 
the SF-121. (The characteristic is indigenous to any VATOL aircraft employing 
aft end thrust deflection for moment control.) X e initially incr eases 
before proceeding to fallow the desired response: i.e. the X g ani., 

likewise, response to pitch attitude are initially in the wrong 
direction. ffle reason for this is that to increase pitch angle the thrust 
must deflect forward (in a di^e-tion to increase X g ) initially to produce a 
nose up pitch mcment. This is a non-minimum phase control characteristic and 

9 

is indicated during control system analysis by zeros of the x t / *LNGSTK 
transfer function lyi'ig in the right half of the s-plane. Thus the SF-121 has 
the potential for i. (Stability and PIO (pilot induced oscillations) if the 
pilot aggressively pursues position control of the airplane. The addition of 
pitch RCS jets acting as a force couple would el ,nate ‘.his tendency and 
allow pitch thrust "<eflection to be used as a direct force control. 
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The longitudinal stick doublet response also demonstrates that there Is a 
modicum of lateral coupling but fairly significant through anticipated heave 
coupling. Thrust deflection reduces the vertical thrust component and the 
aircraft begins to sink (Z e > o). The heave control system provides 
appropriate thrust corrections to drive the heave rate back to trim (Z . 

0). Note that the aircraft will sink for pitch control Inputs In either 
direction from trim. Should It prove annoying the heave coupling can be 
alleviated by a crossfeed from « PITCH to thrust command. 


3.8.5 Case 5 - Lateral Stick Doublet at V A . 10 kt 


The lateral stick doublet response Is the lateral dual of the longitudinal 
stick doublet response. The following verbal and symbology replacements 
inserted in Section 3.8.4 will make the discussion generally applicable as well 
to the time trace of figure 3-28: 


Repl ace 



in Section 3.8.4 with 



The lateral response differs as a dual from the longitudinal response only in 
that Y and subsequently b» do not return to zero (trim) at the end of the 
maneuver whereas X g> the longitudinal dual, does. The difference stems from 
the fact that the fuselage tilts to the right farther and longer than it does 
to the left during the maneuver (i.e. r^ is positive for a longer time 
than it is negative and attains larger positive magnitudes). Since Y e is 
proportional to the Integral of r^ it should, and does, have a positive 
value at the end of the maneuver. Since the aircraft is now moving obliquely 
through the air mass, it also has a steady state 6. 
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Figure 3 - 28 . lateral Stick Doublet Response at V A - IQ kt (Sheet 1 of 3 ) 
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Figure J-28. Literal Stick Doublet Response at V A - 10 kt (Sheet 3 of 3) 


3.6.6 Case 6 - Pedal Pulse at « 10 kt 

The pedal pulse commands the airplane to rotate 28 deg around the Xg 

axis at a rate of approximately 20 deg/sec. This maneuver appears to the 

pilot as a cockpit yaw to the left. As indicated by figure 3-29, the maneuver 

is performed with no control saturation, negligible coupling in the yaw degree 

of freedom, and noticeable, but expected, coupling in the heave, surge, and 

pitch degrees of freedom. The roll control power for this maneuver is 

provided by the RCS jets. The RCS has a demand bleed arrangement and thus, 

when roll control power is required, thrust will decrease. In addition the 

RCS jets are located aft of the eg and thrust along the airplane positive Xg 

axis such that a nose up moment and small positive surge force (negative 

n ) are generated. The reduction in thrust produces a decrease in 
*cg ,• , • 

n x and establishes a sink rate (Z e > 0). The Z e is corrected by the 

hea9e control system. Similarly the pitch moment effects are cancelled by the 

pitch control system. The surge force is not controlled and the airplane 

drifts to a slightly higher airspeed following completion of the maneuver. 

Note that the direction of heave, surge, and pitch coupling is the same 

regardless of the direction of roll control power application. 

In general, the longitudinal coupling displayed by the pedal pulse 
response is typical of RCS - equipped airplanes. Should it prove annoying, 
the culprit can be eliminated or alleviated by one or more of the following 
modifications: 

1. Alleviate pitch coupling by relocating the RCS jets nearer the eg or 
by a s R0LL to «pjjQ R crossfeed. 

2. Eliminate heave coupling by a continuous bleed RCS but at the expense 
of reduced thrust capability. 

3. Alleviate heave coupling by a crossfeed of « R0LL to thrust command. 

3.8.7 Case 7 - Heave Rate Controller Doublet at V A . 10 kt 

The heave rate controller doublet commands the airplane to establish a 10 
ft/sec rate of climb (Z , -10 ft/sec) followed by a 10 ft/sec sink rate 
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Figure 3-29. Pedit Pulse Response at V A • 10 kt (Sheet 3 of 3) 
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(Z fi ■ 10 ft/sec) and return to trim (figure 3-30). The most significant 
feature of this response is the thrust dynamics shown by the trace of the 
thrust before (B4) correction for RCS coupling. Since time constant increases 
and the deceleration limit decreases as thrust level decreases the thrust 
response is slower when, at the same trim speed, establishing a sink rate, 
arresting a climb, or changing from climb to sink than it is when establishing 
a climb, arresting a sink, or changing from sink to climb. 

Coupling into the lateral degrees of freedom was negligible during the 
maneuver and not worth showing in the traces. Coupling into the surge and 
pitch degrees of freedom is small but noticeable; the appropriate trac< j have 
been included. Overall, the response is smooth, stable, and essentially 
single degree of freedom. 

3.8.8 Case 8 - Mode Switching Transient 

To investigate the transients which occur at the switch at 60 kt airspeed 
from conventional t.o hover FCS modes, the airplane was trimmed at 62 kt in a 
0.5 mi radius right turn while decelerating at O.lg along a constant altitude 
flight path. At 0.3 sec into the time history the stick was deflected aft and 
right to impose simultaneous 10 deg/r.ec positive pitch and roll rate commands 
on the airplane. These stick commands were held for 0.7 sec and released and 
should increase airplane pitch and roll angles by approximately 7 degrees 
each. The traces for this maneuver are shown in figure 3-31. 

The 4 and e traces show that the desired maneuver was performed and not 
cancelled or altered by changing FCS modes. Control saturation is displayed 
only by *pj jch* ** exceeds -1.0 for a short time following the release of 
the stick input. The mode switch occurs approximately 1.3 sec into the run. 

It is marked by rapid changes in «rqLL* 4 YAW» 4 PITCH» n y * and 

thrust command (i.e. the heave control system is activat£8) with those in 

4 yaw* n v » and thrust command being fairly large. The pilot will 
. " * 00 . 

likely notice the n^ f which reverses from -.05g to .05g, and «yaw» 
which commands a yaw c 8cceleration step approximately -0.5 rad/sec z to be 
applied to the airplane, since these are unexpected accelerations. He will 
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Figure 3-3Q. Keeve Rite Controller Doublet Response it V A - 10 kt {Sheet I of 2) 
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also notice the thrust Increase to arrest the sink rate established by the 
stick inputs. Piloted simulation studies should determine whether these 
switch transient effects are annoying. For now it sufficies to say that the 
transients are stable, controllable, and, with the few exceptions noted above, 
relatively smooth and predictable. 

3.8.9 Case 9 - Pseudo-Pilot Flown Transition 

The pseudo-pilot flown transition (figure 3-32) demonstrates the capabili- 
ties of the pseudo-pilot transition equations detailed in Section 2.10.2 and 
provide further substantiation of the stability and controllability of the 
airplane through the mode switch region. The transition is initiated by 
retarding the throttle to idle. Psuedo-pi lot's longitudinal stick commands 
are basically open loop in that they are calculated with the aid of a table of 
trim pitch angles vs airspeed. The stick commands are initially proportional 
to pitch rate and become proportional to pitch angle at the mode switch 
point. Psuedo-pi lot's throttle variations are similarly open loop being 
determined from a table of trim throttle settings vs airspeed. Even though 
pseudo-pilot has no inputs to the roll or yaw control systems, the transition 
is well behaved laterally. It can thus be inferred that the baseline FCS will 
enable the transition to be flown with low pilot workload. 

The mode switch, which occurs approximately 56 sec into the transition, is 
accompanied by noticeable rapid changes in « PITCH , « yAH , «R0LL» and 

n y . The change in fipucH is e *P®cted because of the change of command 
fuS<?tion of the longitudinal stick while those in « yAWt and 

n y are unexpected. Thrust command, and consequently, n x , also 
ch£8ge rapidly, not unexpectedly, because the heave con trof 3 sys tern has been 
activated and immediately attempts to arrest the 10 ft/sec sink rate which has 
developed. As discussed in Section 3.8.8 any annoyance related to these 
switch transient effects can best be judged during a piloted simultion. 

3.8.10 Case 10 - Pseudo-Pilot Flown Turn Over a Spot in a 35 kt Wind 

The pseudo-pilot flown turn over a spot in a 35 kt wind (figure 3-33) 
demonstrates the capability of the pseudo-pilot stationkeeping control laws 
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FI pure 3-33. Pseudo-Pilot Flown Tu. n Over i Spot In « 35 kt Wind (Sheet 1 of 5) 
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flgor* 3-33. Pwudo-Pllot Flown Turn Ovtr * Spot In a 35 kt Wind (Shaft 2 of 5) 
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Figure 3-33. Pseudo-Pilot Flown Turn Over e Spot In t 35 kt Wind (Sheet 3 of 5) 
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depicted In figure 2-32. Prior to calculating the time history of the turn, a 
linear analysis was performed of pilot-closed position loops in a headwind to 
establish the pseudo-pilot stationkeeping gains. Because of the non-minimum 
phase characteristics of surge and sway control of the airplane, the gains had 
to be chosen carefully to avoid position loop Instability. The gains 
established by the analysis and used to generate the traces of figure 3-33 are 
as follows: 

K z « -0.006 units *LM5STK^ t; *w * 5.0 ft/ft/sec; 
e s 

* 0.006 units ^v * 5.0 ft/ft/sec j 

e s 

K x - 2.0 ft/sec/ft; * 0. 

The airplane handling characteristics are functions of wind direction but, 
as noted, pseudo-pilot gains were established only for headwind conditions. 

For this reason plus the fact that wind direction changed fairly rapidly 
during the turn it was expected that pseudo-pilot would have a difficult task 
to maintain X g and Y e at the desired X e * Y e * 0 reference point. 

This was indeed the situation. In the time of one complete turn (approxi- 
mately 18 sec), the airplane had drifted in X g f r0 m 10 ft forward to 3 ft 

aft of the reference and was correcting back towards the reference while in 

Y e it had drifted from 10 ft left to 40 ft right of the reference and was 
not yet correcting back towards the reference. Psuedo-pilot had the stick 
properly placed - forward and left - at the completion of the turn for 

correcting both situations. Because there was no wind component along the 

Z e axis, pseudo-pilot did an excellent job, maintaining Z e within 1.5 ft 
of the Z g * -100 ft reference. 

Piloted simulation studies will be beneficial to evaluate the airplane and 
control system performance demonstrated here. At this point it can only be 
concluded that turns over a spot at fairly high turn rates are controllable 
with modest control power usage and have predictable performance. 



The roll rate and « R0LL traces of figure 3-33 demonstrate unexpected 
small one cycle oscillations at 5 and 14 seconds into the maneuver. These are 
postulated to be related to the high « and large e flow conditions which occur 
at these times. Time was not available to study completely this anomaly. 
Fortunately it does not detract from the overall utility of the test case but 
should be studied further to assess its root cause. 

3.8.11 Summary of Results of Performance Evaluation 

The following results are the most significant of the FCS performance 
evaluation: 

1. No stability or controllability problems were actually encountered in 
the test cases. The following potential problems, however, were 
indicated: 

a) Sideslip limitation due to roll and/or pitch control power 
limitations at V ft * 120 kts. (Case 3) 

b) Potential PIO situation due to inherent non-minimum phase 
control characteristics when the pilot attempts tight 
horizontal plane position control in hover. (Cases 4 and 5) 

2. There are noticeable uncommanded mode switching transients 
particularly in the lateral degrees of freedom. (Cases 8 and 9) 

3. There is noticeable coupling into the longitudinal degrees of freedom 
when the RCS is used. (Case 6) 

4. The pseudo-pilot transition and stationkeeping control logic works. 
(Cases 9 and 10) 

5. An unexpected oscillation in roll rate occurs under conditions of 
high o, large s, and relatively high steady state roll rate. Time 
was not available to determine the root cause of this peculiarity. 

The fact that its effect is minor and occurs under seemingly unique 
conditions place it among items recommended for future study. (Case 
10 ) 
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The conclusion drawn fro* the FCS performance evaluation 1$ that the baseline 
FCS provides an adequate starting point for piloted simulation studies and 
that, consequently, the design does not have to be Iterated through the 
performance loop of the FCS design procedure (figure 3-1). 
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4.0 CONCLUSIONS ANO RECONCNDATIONS 


The following conclusions are supported by the developments and analyses 
described herein: 

1. A useful capability has been developed for conducting both piloted 
and non-piloted simulations of the terminal operations of VATOL 
airplanes. This simulation capability can be, but is not restricted 
to being, applied at conceptual design phases where relatively few 
configuration - specific data are available. It will highlight 
handling qualities characteristics which are indigenous to the 
various designs. 

2. The aerodynamics math model is deterministic, functions with OATCOM - 
type data, and can adequately represent the low speed high angle of 
attack, large sideslip aerodynamic characteristics of VATOL airplanes. 

3. The Vought SF-121 airplane, as modeled herein, has adequate control 
power for meeting the trim and maneuvering requirements of the 
MIL-F -83300 and AGARD 577 flying qualities specifications with some 
residual for gust regulation. 

4. The baseline FCS developed for the SF-121 provides an adequate 
starting point for piloted simulation studies. 

5. The most distinguishing flying qualities characteristic of the SF-121 
airplane (ar.d similarly thrust deflected controlled VATOL airplanes) 
is a strong non-minimum phase control characteristic (i.e. initial 
acceleration in wrong direction) which will become apparent when tne 
pilot attempts tight horizontal plane position control in hover. 

The developments and analyses herein have provided background for the 
following recommendations for possible math model modifications and areas of 
emphasis for piloted simulations: 

1. Evaluate the use of rotated inertial axis system coupled with Euler 
angles to orient the airplane body axes (Section 2.9). This 
alternative to the direction cosine formulation will circumvent the 
singularity in the standard Euler transformation at e « 90 deg. and 
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provide continuous attitude references for attitude control loops 
throughout the operating range of VATOL airplanes. If the approach 
has merit, incorporate the appropriate relations into the math modH. 

2. Determine and correct, if necessary, the mechanism for the small 
oscillation in roll rate which appears to occur under conditions of 
simultaneous high angle of attack rates, large sideslips, and 
relatively high steady state roll rates (Section 3.8.10). 

3. Piloted simulation studies should consider at least the following 
issues regarding terminal nper»tions of VATOL airplanes: 

o Non minimum phase control characteristics and their impact 
on pilot comfort, workload, and ability to effect precise 
position control in hover and on the need for additional 
moment controls (e.g. independent pitch and yaw RCS). 
o Cockpit and/or pilot rotation during transition. First 

consider its necessity. If required, then consider whether 
switching of cockpit controller roles and pilot/cockpit 
angle control should be automatic or manual functions, 
o Transition control system considerations - control system 
type (rate command, rate command/attitude hold, attitude 
command, etc), automatic or manual mode switching, mode 
blending, annoyance level of uncommanded motions during 
mode switching, co^kit controller authorities and 
sensitivities, pilot input command shaping requirements, 
system bandwidths, gust regulation capabilities. 
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